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Chapter  1 
Introduction 


In  recent  years,  new  complex  structures  have  been  fabricated  on  the  nano-scale  which 
exhibit  novel  optical  properties.  An  example  of  such  devices  include  Photonic  crystals 
(PhC).  A  PhC  is  a  periodic  refractive  index  perturbation  in  a  high-index  material.  These 
man-made  structures  can  forbid  light  propagation  within  certain  frequency  ranges  and 
demonstrate  novel  dispersion  characteristics  (1). 

In  order  to  fabricate  such  complex  structures,  new  techniques  have  been  developed; 
one  particular  method  involves  particle  self-assembly.  In  the  self-assembly  method  (2) 
particles  are  deposited  upon  a  substrate  and  organize  themselves  in  a  form  determined 
by  their  inter-particle  interactions.  A  fundamental  constraint,  however,  in  the  self- 
assembly  method  is  that  only  energetically  favorable  structures  are  formed.  Templated 
self-assembly  methods  have  been  introduced  (3;  4)  to  expand  the  range  of  structures 
that  can  be  fabricated  using  self-assembly.  In  this  method,  particle  assembly  occurs 
on  a  template.  The  geometry  of  the  template  guides  particle  arrangement,  allowing 
the  formation  of  diverse  crystal  structures.  A  fundamental  problem  with  this  method 
is  that  the  final  structure  is  determined  strongly  by  the  geometry  of  the  template. 
Consequently,  the  formation  of  structures  that  possess  different  size  particles  and  a 
wide  range  of  material,  is  not  well  controlled. 

On  the  other  hand,  since  the  discovery  of  optical  tweezers  (5;  6),  optical  forces  near 
micro-photonic  devices  have  been  studied  extensively.  While  standard  optical  tweezers 
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use  a  focused  laser  beam  to  trap  and  manipulate  particles,  structured  light  fields  allow 
even  greater  control.  Large  optical  forces  have  been  achieved  by  using  the  strong  elec¬ 
tromagnetic  fields  near  nano-apertures  (7),  gratings  (8),  photonic-crystal  micro-cavities 
(9;  10),  slot  waveguides  (11)  and  plasmonic  structures  (12;  13).  Applications  of  such 
particle  manipulation  have  been  found  from  the  physical  sciences  to  biology  (14;  15). 
However,  much  of  this  work  has  only  considered  single-particle  traps. 

In  this  dissertation,  I  study  optical  forces  near  micro-fabricated  devices  for  multi¬ 
particle  manipulation.  I  consider  particles  of  different  sizes  and  compositions.  In 
particular,  I  focus  my  study  on  both  dielectric  and  gold  particles  as  well  as  Giant 
Unilamellar  Vesicles  (GUVs). 

First,  I  consider  optical  forces  near  a  PhC  and  establish  the  feasibility  of  a  technique 
which  we  term  Light-Assisted  Templated  Self-assembly  (LATS).  In  contrast  to  previous 
work  on  Fabry-Perot  enhancement  of  trapping  forces  above  a  flat  substrate  (16),  I 
exploit  the  guided  resonance  modes  of  a  PhC  (17)  to  provide  resonant  enhancement  of 
optical  forces.  The  guided  mode  forms  spatially  distinct  trapping  patterns  at  different 
polarizations  and  wavelengths.  I  consider  the  possibility  of  assembling  multi-particle 
patterns.  Specially,  I  show  how  optical  forces  near  a  photonic  crystal  slab  can  be  used 
to  assist  templated  self-assembly  and  permit  different  particle  pattern  formation.  Thus, 
a  PhC  will  serve  as  an  optically  reconfigurable  template  for  multi-particle  trapping. 

Furthermore,  in  contrast  to  localized  mode  trapping  techniques  (which  is  necessarily 
limited  to  a  smaller  area,  and  thus,  fewer  trap  sites)  the  use  of  an  extended  mode  has 
obvious  advantages.  The  fact  that  every  trap  operates  simultaneously  in  a  small  area 
eliminates  the  need  for  a  spatial  light  modulator  or  beam  scanning  of  the  sample.  This 
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is  a  unique  feature  of  our  extended-mode  trapping  system. 

Combining  the  optical  field  created  in  our  system  with  a  microfluidic  flowalso 
provides  the  potential  for  sorting  particles  via  the  creation  of  lock-in  states  (18). 
Additionally,  changing  the  shape  of  the  hole  within  each  unit  cell,  an  asymmetric 
potential  could  be  created  which  could  operate  as  a  ratchet  (19).  We  anticipate  that 
our  system  could  be  applied  in  self-assembly  experiments  and  also  to  build  clusters  in 
a  reversible  way;  the  later  would  enable  the  creation  of  reconfigurable,  optical  matter 
(20).  Our  LATS  system  is  able  to  assemble  structures  with  arbitrary  symmetries  which 
are  ordinarily  unavailable  for  assembly  due  to  energy  constraints.  We  demonstrate  this 
unique  ability  both  theoretically  and  experientally.  In  addition,  our  structure  could  be 
fabricated  upon  an  active  material,  thereby  merging  the  source  and  trap.  With  proper 
electrical  excitation,  such  structures  could  easily  be  adopted  for  various  applications. 
Finally,  although  previous  work  has  examined  light  forces  between  different  layers  of 
particles  in  colloidal  PhC  (21;  22),  forces  on  colloidal  layers  above  a  micro-fabricated 
template  were  not  considered.  I  used  full-vectorial  electromagnetic  simulations  of 
Maxwell’s  equations  in  time  domain  to  calculate  optical  forces  on  dielectric  particles  of 
500  nm. 

In  Chapter  3,  I  present  the  experimental  demonstration  of  LATS.  I  explain  all 
relevant  setup  and  fabrication  methods.  Next,  I  use  a  particle  tracking  software  to 
record  the  position  of  a  520  nm  polystyrene  particle  while  trapped  in  a  square  lattice 
PhC.  In  addition,  I  numerically  calculate  the  optical  force  at  the  center  of  the  trapping 
pattern  and  the  corresponding  optical  stiffness  of  the  PhC  trap. 
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In  the  following  chapter,  I  calculate  the  optical  force  near  a  PhC  micro-cavity 
for  gold  and  dielectric  particles  ~  50  nm  in  diameter  in  diameter.  Next,  I  use  the 
electromagnetic  field  re-distribution  around  a  trapped  gold  nano-particle  to  direct  a 
metal-dielectric  cluster  formation. 

In  the  final  chapter,  I  explore  optical  forces  near  a  Dual  Beam  Optical  Trap 
(DBOT)  and  measure  the  bending  modulus  of  a  GUV  ~  10  /im.  First,  I  present  a 
method  to  extract  the  bending  modulus  of  the  membrane  from  the  area  strain  data. 
This  method  incorporates  three-dimensional  ray-tracing  to  calculate  the  applied  stress 
in  the  DBOT  within  the  Ray  Optics  approximation.  I  compare  the  optical  force 
calculated  using  the  Ray  Optics  approximation  and  Maxwell  Stress  Tensor  method  the 
ensure  the  approximation’s  accuracy.  Next,  we  apply  this  method  to  3  populations 
of  GUVs  to  extract  the  bending  modulus  of  membranes  comprised  of  saturated  and 
monounsaturated  lipids  in  both  gel  and  liquid  phases.  We  find  that  bending  modulus 
values  depend  on  lipid  phase  and  are  normally  distributed  for  the  three  populations  we 
tested.  Our  results  also  indicate  that  the  addition  of  cholesterol  to  POPC  membranes 
does  not  alter  the  respective  bending  modulus.  Moreover,  the  values  we  obtain  for 
bending  moduli  agree  well  with  those  found  in  literature.  This  technique  is  a  promising 
route  to  obtaining  statistically  relevant  bending  moduli  data  for  lipid  membranes  with 
controlled  compositions. 


1.1  Forces 

For  a  particle  in  water,  we  can  consider  the  interplay  of  various  forces:  gravitational, 
buoyant,  brownian,  electrostatic  and  those  forces  that  develop  from  the  interaction 
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between  an  external  electromagnetic  field  and  a  particle.  I’ll  refer  to  the  last  one  as  an 
optical  force. 

In  general,  the  dynamics  of  the  particle  is  governed  by  the  following  equation: 


UIX  'yX  Ffluid  “f  Foptical  F  X(f)  (1-1) 

where  7  is  the  friction  coefficient,  and  are  mutually  uncorrelated  white  noises 
that  obey  the  fluctuation-dissipation  relation  (Xi(f)IXi(f))  =  2kBT'ySijS{t  —  t'). 

In  this  chapter,  I  describe  the  origin  and  order  of  magnitude  for  these  forces  acting 
upon  a  I  /im  polystyrene  particle. 

1.1.1  Mechanical  forces 

In  the  particular  case  of  a  polystyrene  particle  in  water  {plater  =  ^9/cm^,ri^ater  = 
1.3  X  Ppoiy  =  2.5g/cm^),  the  gravitational  force  is  on  the  order  of  1  x  10“^ 

pN  while  the  buoyancy  force  is  ~  0.6  x  10“^  pN.  Thus,  the  particle  feels  a  net 
gravitational  force  on  the  order  of~  1  x  10“^  pN. 

Using  the  diffusion  coefficient  D  =  ksT/Qiirri  =  0.5/xm^/s  as  well  as  the  drag 
coefficient  7  =  Qnrpr],  the  average  Brownian  force  can  be  shown  to  be  ~  1  X  10  ^  pN. 

For  this  system  the  Reynolds  number  is  ~  1  x  10“^.  At  a  low  Reynolds  number, 
the  inertia  plays  no  role  in  the  particle  dynamics  whatsoever.  That  means,  what  the 
particle  is  doing  at  any  given  moment  is  entirely  determined  by  the  forces  that  are 


5 


17 


exerted  on  particle  at  that  moment,  and  by  nothing  in  the  past.  Thus,  Eq.  (1.1)  reduces  to 


yX  T/Zmd  T  Foptical  T  X.it') 


(1.2) 


1.1.2  Van  der  Waals  force 

The  Van  der  waals  force  arises  from  electron  cloud  fluctuations  surrounding  the  nucleus 
of  an  electrically  neutral  atom.  This  force  is  a  short  range  force  and,  in  generally,  is 
attractive. 


For  macroscopic  bodies  with  known  volumes  and  numbers  of  atoms  (molecules 
per  unit  volume),  the  total  van  der  Waals  force  is  often  computed  based  on  the  micro¬ 
scopic  theory  as  the  sum  over  all  interacting  pairs.  It  is  necessary  to  integrate  over  the 
total  volume  of  the  object,  which  makes  the  calculation  dependent  on  the  object’s  shape. 

The  van  der  Waals  force  between  two  spheres  of  constant  radii  Ri  and  i?2  is  then  a 
function  of  separation  r 

^  AR1R2 

=  ~erHR,  +  R,) 

where  A  is  a  constant  that  depends  on  both  materials.  For  two  500  nm  polystyrene 
particles  in  water,  A  =  0.22  x  10“^°  J  and  F  =  1  x  10“^  pN  at  300  nm  of  separation. 


In  the  case  of  a  sphere  near  a  surface,  the  van  der  Waals  force  is 

_AR 

^VW  — 


(1.4) 
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For  a  500  nm  particle  floating  in  water  100  nm  from  a  silicon  surface  F  =  1  x  10  ^  pN. 

1.1.3  Optical  forces 

Optical  forces  arise  whenever  an  electromagnetic  wave  propagates  in  a  medium.  It 
can  be  understood  as  a  change  in  momentum  of  the  propagating  light  due  to  a  change 
in  refractive  index;  it  is  a  force  felt  by  the  induces  dipoles  in  the  medium  due  to  the 
presence  of  an  electromagnetic  field. 

Maxwell  stress  tensor  method 

For  a  small  volume  V  the  induced  charge  p  will  experience  a  mechanical  force  given  by 
the  Lorent’s  relation  (23): 


F=  {pE  +  JxB)dV  (1.5) 

Jv 

Expressing  the  charge,  current  and  magnetic  field  in  terms  of  the  total  electric  field  E 
we  can  obtain 

F=  I  f^-dd-^  (1.6) 

Jss  dt 

where  55'  is  the  surface  enclosing  the  volume  V,  S  is  the  Poynting  vector  and  Tap  is  the 
Maxwell  stress  tensor  defined  as 

Tap  =  DaEp  —  -5ap{D  ■  E)  +  HaBp  —  -5ap{H  ■  B)  (1-7) 

In  general,  there  are  two  natural  approximations  often  employed  in  such  calculations: 
particles  much  smaller  than  the  wavelength  (dipole  approximation)  or  much  larger  (ray 
optics). 
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Optical  force  using  dipole  approximation 

Let  us  write  Eq.  (1.6)  in  the  volume  form  (24;  25) 

{F)  =  I  {f^)-da=  [  dVV(T(r))  (1.8) 

Jss  J 

For  a  small  particle,  the  total  field  inside  the  particle  can  be  expressed  in  terms  of  the 
external  field  Eo{r}  as 

E{r)  =  Eo{r)  +  aG{r  -  ro)Eo{r)  (1.9) 


(1.10) 


(1.11) 


The  time  average  force,  expresed  in  terms  of  the  external  electric  field  and  the  induced 
dipole  moment  is  (26) 

d.EoiVE'oj  (1.12) 


{F)  =  -Re  J2p^VE*,  =-Re 
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using  Faraday’s  law  V  x  E  =  iu^oH  and  the  relation  EiVE*  =  {E  ■  V)E*  +  E  x 
(V  X  E),  the  dipolar  force  (F)  can  be  expand  as  the  sum  of  the  gradient  force,  radiation 
pressure  and  spin  density  of  the  light  field  as 


(F) 


1 

2 


Re  [aEo  ■  V-Eq]  +  o' scat 


+  O-scaticV  X  {Ls)[ 


(1.13) 


where  ascat  =  koIm{a)/eo  is  the  scattering  cross  section  and  (Ls)  =  ^(Fq  x  Eq)  is 
the  time  averaged  spin  density. 


The  first  term  is  the  gradient  force,  the  second  term  is  the  radiation  pressure  (pro¬ 
portional  to  the  Poynting  vector  (S)),  and  the  last  term  is  a  curl  force  associated  to  the 
nonuniform  distribution  of  the  spin  density  of  the  light  field. 

Optical  binding 

The  dipole  approximation  done  in  Eq.  (1.13)  is  valid  for  a  single  particle  in  an  external 
field.  For  multiple  particles,  the  total  electric  field  is  the  superposition  of  the  incident 
electric  field  and  the  field  emitted  by  the  induced  dipoles. 

Let  us  consider  two  dielectric  particles  with  induced  dipole  moments  p(rA)  and 
p(rs)  in  an  electric  field  E  such  that,  Pi(rAiB)  =  ct^-^Ej(rAiB)-  In  terms  of  the  incident 
field  Eq,  the  total  field  can  be  rewritten  as  (27): 


Ei{rA)  =  Eoi{rA)  +  Gij{rA,  rB)af^Ek{rB)  (1-14) 

Ei{rB)  =  EQi{rB)  +  Gij{rB,rA)afkEk{rA)- 

The  solution  for  particle  B  of  this  equation  is  Ei{rB)  =  KijE^jirB)  + 
KijGjka^^EQrni^A),  where  Kij  is  the  inverse  tensor,  [I  —  G{rB-,  r A)oi^G{r A-,r B)oi^]~^ ■ 
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For  small  particles  we  can  approximate  Kij  =  5ij,  thus,  Eq.  (1.12)  for  particle  B  can  be 
written  as 

F.irs)  =  ]^Re[o^*El^{rB)d^Eo^{rB)  (1.15) 

{rA)d^Eoi  {tb  ) 
+<^ij*F*oj{rB)d^Giqaq^Eon{rA)]- 


The  first  term  is  the  classical  formula  for  the  optical  force  which  is  determined  solely 
by  the  spatial  distribution  of  the  incident  field.  This  force  is  usually  further  divided 
into  both  a  gradient  force  and  scattering  force.  The  second  term  is  the  interaction 
between  the  incident  field  and  the  induced  dipole  from  the  scattered  field.  The  third 
term  corresponds  to  the  force  acting  upon  the  induced  dipole  from  the  scattered  field  by 
the  incident  field. 

In  general,  by  using  the  Maxwell  stress  tensor  method  we  incorporate  the  effect  of 
optical  binding. 

Optical  force  in  the  ray  optics  approximation 

For  particles  much  larger  than  the  wavelength  of  the  incident  light,  we  can  use  the 
Ray  Optics  approximation  (RO).  Under  this  assumption,  rays  undergo  a  change  in 
momentum  after  a  reflection  on  a  surface.  This  change  in  momentum  induces  an  optical 
force  on  the  object. 

Let  us  consider  a  ray  propagating  along  the  incident  direction  iq.  This  ray  carries  a 
momentum  of  p  =  {ni/c)PiQ,  where  ni  is  the  refractive  index  of  the  medium,  P  is  the 
optical  power  and  c  is  the  velocity  of  light  in  vacuum.  Following  a  reflection,  the  force 
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acting  upon  a  surface  can  be  calculated  as  the  difference  in  momentum  carried  by  the 
incident  ray  compared  to  that  carried  by  both  the  reflected  and  transmitted  rays  (28;  29), 


/ 


niP 

c 


*0  ~  {Ro'f' 0  H - Toio) 

rii 


(1.16) 


where  foand  to  define  the  unit  vectors  along  the  reflected  and  transmitted  directions, 
respectively.  Rq  and  To  represent  the  Fresnel  reflection  and  transmission  coefficients. 


For  optical  fibers,  the  intensity  of  light  emanating  is  best  approximated  by  a  Gaussian 
distribution.  Thus,  the  power  incident  on  a  differential  area  As  is  related  to  the  total 
beam  power  by  the  following  equation. 


P 


As  cos6 


(1.17) 


where  w  and  r  represent  the  beam  width  and  the  radial  distance  from  the  beam  axis  i 
at  the  interacting  point,  and  9  is  the  angle  between  the  surface  normal  and  the  beam 
propagation  direction.  Using  Eq.  (1.16)  and  (1.17)  the  surface  stress  can  be  shown  to  be 


a  = 


riiCosO 


*0  ~  {Rq^ 0  H - Toio) 

ni 


2Pt  (  2P 

— -exp - ^ 

TTW^  V 


(1.18) 


In  the  case  of  a  close  object,  a  transmitted  beam  experiences  multiple  internal  reflections. 
If  we  consider  m  internal  reflections,  the  total  optical  surface  stress  is 


^  2PTniToAscos6 

o-m  =  - - exp 

nw^cAsm 


2P 

w‘^ 


m—1 

i=i 


_  (—R  r  +T  i  ) 

V  m  '  m^m ) 

rii  Til 


(1.19) 


The  incident  vector  takes  into  account  the  divergence  of  the  beam,  meanwhile  the 
subsequent  reflected  and  transmitted  vectors  can  be  computed  using  Snell’s  law.  Finally, 
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the  total  force  is  the  integral  of  the  stress  over  the  surface. 


1.1.4  Thermal  Gradient  force 

In  a  material  that  possesses  a  complex  dielectric  constant,  the  absorption  of  electromag¬ 
netic  energy  induces  a  change  in  temperature  which  leads  to  a  convective  force  due  to 
the  fluid  flow. 

The  dynamics  of  such  a  system  can  be  studied  by  numerically  solving  the  steady- 
state,  incompressible  Navier-Stokes  equation  in  the  microfiuidic  environment  together 
with  a  heat  transfer  equation.  The  incompressible  Navier-Stokes  equations  consist  of  a 
momentum  balance 


p—  +  pu  ■  Vm  =  — Vp  -f  riV'^u  +  F  (1-20) 

ot 

This  equation  ignores  variations  in  density  with  temperature.  We  add  free  convection 
to  the  fluid  flow  with  the  Boussinesq  approximation:  F  =  —g{p{T)  —  p(To)).  We 
assume  that  the  fluid  density  p  is  independent  of  pressure  and  depends  linearly  on  the 
temperature  T:  p  —  po  =  ~PoC'p(T  —  Tq),  where  Cp  is  the  heat  capacity. 

The  heat  equation  is  an  energy  conservation  equation.  The  change  in  energy  is  equal 
to  the  heat  source  minus  the  divergence  of  the  diffusive  heat  flux 


pCp 


+  U-VT 


+  V  •  i-kVT)  =  Q 


(1.21) 
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where  Q  is  a  source  term.  The  solution  of  the  coupled  system  provides  the  velocity 
field  and  temperature  distribution  within  the  fluid. 

This  effect  is  important  for  systems  with  absorption.  A  numerical  simulation  esti¬ 
mates  that,  for  a  1  /i  m  particle,  the  convective  force  is  ~  1  x  10“®  pN  for  the  channel 
and  parameters  used  in  this  thesis. 
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Chapter  2 

Light  Assisted  template  self-assemhly 


Self-assembly  methods  (2)  have  been  used  to  eonstruet  eomplex  materials  ineluding 
three  dimensional  photonie  erystals  (30).  However,  a  fundamental  eonstraint  is  that  only 
energetieally  favorable  struetures  are  formed.  To  overeome  this  problem,  templated 
self-assembly  methods  have  been  introdueed  (3;  4).  The  template  guides  partiele 
arrangement,  allowing  the  formation  of  diverse  albeit  statie  erystal  struetures.  Here 
we  eonsider  how  optieal  forees  ean  be  used  to  direet  assembly  and  reeonfiguration 
of  partieles  on  a  photonie  erystal  slab,  whieh  serves  as  an  optieally  reeonfigurable 
template.  Our  ealeulations  prediet  the  formation  of  stably  trapped  erystal  patterns 
that  depend  on  the  wavelength  and  polarization  of  the  light  souree.  We  envision  that 
the  proeess  of  light- assisted  templated  self-assembly  may  be  used  for  fabrieation  of 
eomplex  photonie  materials  and  all-optieally  reeonfigurable  photonie  deviees. 

Optieal  forees  have  been  studied  extensively  sinee  the  diseovery  of  optieal  tweezers 
(5;  6).  While  standard  optieal  tweezers  use  a  foeused  laser  beam  to  trap  and  manipulate 
partieles,  struetured  light  fields  (14;  15)  allow  even  greater  eontrol  over  partiele 
defleetion,  transport,  and  sorting.  Reeently,  researehers  have  used  light  fields  near 
mierophotonie  deviees  to  trap  and  manipulate  partieles.  This  work  leverages  the 
strong  eleetromagnetic  gradients  near  deviees  sueh  as  nanoapertures  (7),  gratings  (8), 
photonie-erystal  mieroeavities  (9;  10)  to  generate  large  optieal  forees.  However,  mueh 
of  this  work  has  eonsidered  single-partiele  traps. 
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Figure  2.1:  Light-assisted  templated  self-assembly  using  a  photonic  crystal  slab 

Here  we  suggest  the  possibility  of  assembling  multi-partiele  patterns.  We  propose 
to  use  the  guided  resonanee  modes  of  photonie  erystal  slabs  (17)  to  enhanee  the  optieal 
forees  aeting  upon  partieles  in  a  solution.  In  eontrast  to  previous  work  on  Fabry-Perot 
enhaneement  of  trapping  forees  above  a  flat  substrate  (16),  the  photonie  erystal  is  used 
to  form  spatially  distinet  trapping  patterns  at  different  polarizations  and  wavelengths. 
While  previous  work  has  examined  light  forees  between  different  layers  of  partieles 
in  eolloidal  photonie  erystals  (22),  forees  on  eolloidal  layers  above  a  mierofabrieated 
template  were  not  eonsidered. 

The  system  we  eonsider  is  shown  sohematieally  in  Fig  2.1.  In  an  experiment, 
light  would  illuminate  the  photonie  erystal  slab  from  below.  The  struetured  light 
fields  above  the  slab  give  rise  to  optieal  forees  on  partieles  in  a  solution,  whieh  ean 
potentially  result  in  trapping.  Due  to  the  variation  of  the  eleetromagnetie  field  over 
spatial  length  scales  comparable  to  the  particle  diameter,  the  dipole  approximation 
cannot  be  assumed.  However,  the  optical  force  can  be  found  by  using  full-vectorial 
calculations  of  Maxwells  equations  to  calculate  the  Maxwell  Stress  Tensor  (22). 
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2.1  Design 

For  concreteness,  we  consider  a  silicon  photonic  crystal  (e  =  11.9)  with  a  square  lattice 
of  holes  of  radius  0.2a  and  thickness  0.6a  resting  on  a  silica  substrate  (e  =  2.1),  where 
a  is  the  lattice  constant  of  the  photonic  crystal.  The  holes  in  the  photonic  crystal  and 
the  region  above  the  crystal  are  filled  with  water  (e  =  1.77).  The  incident  light  is  an 
x-polarized  plane  wave  propagating  in  the  2:  direction. 

We  calculate  the  transmission  through  the  slab  numerically  using  the  finite- 
difference  time-domain  (FDTD)  method  (31)  with  the  freely  available  MIT  MEEP 
package  (32).  The  computational  resolution  is  32  grid  points  per  lattice  constant. 
Eig  2.2(a)  shows  the  normalized  transmission  through  the  photonic-crystal  slab.  The 
transmission  exhibits  three  distinctive  resonance  peaks  at  frequencies  z/(c/a)  =  0.4721, 
0.4958,  and  0.5008.  The  characteristic  shape  of  the  resonance  peaks  is  well  known 
from  existing  work  in  literature  (17).  We  will  refer  to  the  three  resonances  shown  as 
Rl,  R2,  and  R3  respectively. 


2.2  Optical  force  using  guided  modes 

We  calculate  the  force  on  polystyrene  spheres  (e  =  2.28)  of  radius  0.1a.  The  center 
of  the  sphere  is  placed  at  a  height  of  0.25a  above  the  slab  surface.  We  calculate  the 
Maxwell  stress  tensor  on  the  surface  of  an  integration  box  of  size  0.3a  centered  on  the 
particle  (23).  The  calculation  is  done  in  a  computational  cell  that  includes  one  period  of 
the  photonic  crystal  slab.  Periodic  boundary  conditions  are  used  in  the  lateral  directions 
(x  and  y),  and  perfectly  matched  layer  boundaries  are  used  in  the  z  direction.  The  forces 
calculated  in  this  manner  represent  the  optical  forces  on  a  periodic  array  of  polystyrene 
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(a) 


Figure  2.2:  (a)  Normalized  transmission  of  a  silicon  slab  on  a  silica  substrate  for  a  height  of 
0.6a  and  hole  radius  of  0.2a,  (b)  Force  in  the  z  direction  for  a  sphere  of  radius  0.1a  placed  0.25a 
above  the  slab  surface.  The  sphere  is  positioned  equidistant  from  the  four  closest  holes. 

spheres  above  the  photonie  erystal  slab. 

We  plot  the  dimensionless  quantity  F^c/cf),  where  Fz  is  the  z  eomponent  of  the 
optical  force,  (j)  is  the  incident  power  per  unit  cell  and  c  is  the  speed  of  light  (see  Fig 
2.2(b))  for  a  sphere  positioned  equidistant  from  the  four  closest  holes.  For  each  of 
the  three  resonance  frequencies  shown  in  Fig  2.2(a),  we  observe  a  peak  in  the  optical 
force.  The  sign  of  the  force  is  negative,  indicating  that  the  particles  are  attracted  toward 
the  slab.  For  comparison,  we  calculate  the  radiation  pressure  of  a  plane  wave  on  a 
spherical  particle  in  the  absence  of  the  photonic  crystal.  We  obtain  a  repulsive  force  of 
Fzc/(l)  ~  1  X  10  three  orders  of  magnitude  smaller  than  the  attractive  optical  force 
above  the  photonic  crystal.  We  expect  that  even  larger  optical  forces  can  be  achieved  by 
using  a  photonic-crystal  slab  with  higher-Q  resonances.  Higher  Q  has  been  correlated 
with  decreasing  hole  size  (17). 
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Figure  2.3:  Normalized  intensity  for  normally  incident  x-polarized  light  for  Rl,  R2,  and  R3, 
respectively  at  a  height  of  (a-c)  0.25a  above  the  surface  of  the  slab  and  (d-f)  0.1a  above  the 
surface  of  the  slab.  White  circles  indicate  hole  positions. 


The  intensity  for  these  resonanees  is  shown  in  Fig.  2.3(a-e)  at  a  height  of  0.25a 
above  the  surfaee  of  the  slab.  Fig  2.3(d-f)  shows  the  intensity  at  a  height  of  0.1a  (equal 
to  the  partiele  radius).  The  intensity  \E\'^  is  normalized  to  the  souree  intensity  \Eq\^. 
We  observe  that  eaeh  resonanee  has  a  different  field  profile.  There  is  a  slight  ehange  in 
field  profile  with  height,  and  the  intensity  inereases  as  you  approaeh  the  slab. 


2.3  Optical  force  maps  and  reconfiguration 

We  now  look  at  the  spatial  dependenee  of  the  foree  on  resonanee.  The  foree  at  eaeh 
point  is  ealeulated  by  plaeing  the  sphere  at  that  point,  ealeulating  the  full-veetorial 
eleetromagnetie  fields,  and  eomputing  the  Maxwell  Stress  Tensor.  For  numerieal 
reasons,  we  have  ealeulated  the  forees  at  a  height  of  0.25a  above  the  slab  to  insure  that 
the  Maxwell  Stress  Tensor  integration  surfaee  does  not  overlap  with  either  the  partiele 
or  the  slab  boundaries.  Figs  2.4(a),  (b),  and  (e)  show  the  foree  on  a  vertieal  (x-z)  shoe 
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Figure  2.4:  (a-c)  Force  in  the  x-z  plane  at  y  =  0,  (d-f)  Force  in  the  x-y  plane  at  z  =  0.25a  for 
Rl,  R2,  and  R3,  respectively.  The  red  dots  represent  stable  trapping  points  and  the  black  circles 
indicate  hole  positions. 


at  y=0.  The  length  of  the  arrows  next  to  the  graphs  represents  a  value  Fc/<p  =  0.2.  The 
vertieal  forees  are  attraetive  over  the  whole  unit  eell:  this  indieates  that  a  partiele  near 
the  slab  will  be  attraeted  to  it. 


Figs  2.4(d),  (e),  and  (f)  show  the  foree  on  a  horizontal  (x-y)  sliee  0.25a  above 
the  slab  surfaee.  The  eireles  indieate  hole  positions.  The  arrows  next  to  the  graphs 
represent  a  value  of  Fc/0  =  0.2,  and  the  red  eireles  represent  the  stable  points  to  whieh 
the  partieles  are  attraeted.  Rl,  R2  and  R3  eorrespond  to  different  stable  patterns.  Fig 
2.4(d)  has  an  additional  stable  point  at  the  eenter  of  the  unit  eell,  but  it  is  not  labeled 
with  a  red  dot  beeause  it  has  a  weak  lateral  foree  eompared  to  other  trapping  points. 
Similarly,  Fig  2.4(e)  has  two  additional  weakly  stable  points  at  the  bottom  (or  top) 
edge  of  the  unit  eell.  We  note  that  the  spatial  profile  of  the  forees  exerted  on  partieles 
touehing  the  slab  surfaee  (height  0.1a)  may  differ  slightly  from  those  shown  here  due 
to  the  weak  variation  in  intensity  with  height  (Fig  2.3). 
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Figure  2.5:  (a)  Pattern  for  x-polarization,  (b)  Pattern  for  y-polarization. 


Different  light  polarizations  eorrespond  to  different  self-assembled  patterns,  as 
shown  in  Fig  2.5.  Due  to  the  symmetry  of  the  photonie  erystal,  the  patterns  for 
y-polarized  light  are  the  same  as  those  for  x-polarized  light,  only  rotated  7r/2  (Figs 
2.5(a)  and  2.5(b)). 


In  an  experiment,  partieles  may  be  trapped  one  by  one,  or  a  few  at  a  time,  in  order 
to  build  up  ordered  patterns.  We  believe  that  the  red  dots  shown  in  Figs  2.4  and  2.5  are 
good  predietors  of  the  patterns  that  will  be  formed  in  this  manner.  In  our  eomputations, 
due  to  the  boundary  eonditions  on  the  unit  eell,  we  are  ealeulating  the  optieal  foree  on 
periodie  arrays  of  partieles  above  the  photonie  erystal  slab.  However,  sinee  interpartiele 
interaetions  (27)  are  negligible  in  this  system,  the  foree  on  an  individual  partiele  above 
the  slab  will  be  nearly  identieal. 


2.4  Inter-particle  interaction 

To  show  that  interpartiele  interaetions  are  weak  eompared  to  the  foree  on  a  partiele 
from  the  photonie  erystal,  we  eonsidered  the  geometry  shown  in  Fig  2.6(a).  A  first 
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Figure  2.6:  (a)  Schematic  of  the  system  (b)  Transmission  spectrum  of  photonic-crystal  slab  with 
and  without  particles  above  (c)  Fxc/(f>  on  particle  1  for  3  different  positions  of  particle  2  (d) 
Fyc/cj)  on  particle  1  for  3  different  positions  of  particle  2.  The  legend  indicates  the  position  of 
particle  2  in  (b),  (c),  and  (d). 


particle  (labeled  1,  blue  circle)  was  placed  over  the  slab  at  a  position  of  [0.3a,  0.3a, 
0.25a].  A  second  particle  was  placed  in  one  of  three  alternate  positions  (labeled  2, 
green  circles).  Fig  2.6(b)  shows  that  the  presence  of  the  particles  does  not  visibly  affect 
the  transmission  spectrum.  In  Figs  2.6(c)  and  (d),  we  plot  the  lateral  forces  on  particle 
1  for  different  positions  of  particle  2.  It  can  be  seen  that  particle  2  has  minimal  effect 
on  the  force. 


As  an  additional  check,  we  recalculated  the  force  maps  of  Figs  2.4(d-f)  using  a 
larger  computational  cell  of  size  2a  x  2a,  with  one  particle  per  computational  cell.  No 
changes  were  observed  in  the  force  maps. 

Lastly,  we  checked  that  when  particles  were  placed  at  each  of  the  red  dots  in  Fig 
2.4(d-f),  the  configuration  remained  stable  despite  the  perturbation.  One  particle  per 
unit  cell  was  displaced,  and  We  verified  that  the  restoring  force  was  in  the  necessary 
direction  to  restore  the  initial  pattern. 


21 


33 


In  the  system  studied  here,  the  particles  do  not  significantly  affect  the  transmission 
spectrum  (Fig  2.6(b)).  We  have  checked  that  even  if  the  particles  are  placed  inside  the 
holes  in  the  slab  or  on  the  slab  surface,  the  shift  in  the  resonance  is  a  small  fraction  of 
the  resonance  width  (  15%  or  less).  For  larger  particles,  particles  of  higher  index,  and/or 
higher  Q  slab  resonances,  particle  trapping  may  begin  to  shift  the  resonance  position. 
This  opens  up  the  possibility  for  intriguing  phenomena  such  as  self-induced  or  bistable 
trapping,  effects  which  have  previously  been  studied  for  particle  trapping  near  photonic 
cavities  (9).  This  is  an  interesting  direction  for  further  research. 


2.5  Parametric  dependency 

Fig  2.7(a)  shows  the  vertical  force  as  a  function  of  the  dielectric  permittivity  of  the 
sphere,  e,  for  a  sphere  of  radius  0.15a.  The  magnitude  of  the  force  increases  with  e,  and 
no  shift  is  observed  in  the  resonance  peak.  Similarly,  Fig  2.7(b)  shows  the  vertical  force 
as  a  function  of  sphere  radius  for  a  particle  with  e  =  3.  The  force  increases  with  increas¬ 
ing  sphere  radius  and  no  shift  is  observed  in  the  resonance  peak.  For  larger  sphere 
radii,  increasing  e  will  increase  the  force  and  shift  the  resonance  peak,  as  shown  in  Fig 
2.7(c).  The  figure  shows  the  region  of  the  spectrum  near  R1  for  a  sphere  radius  of  0.35a. 

Changing  the  template  will  change  the  frequencies  and  quality  factors  of  the 
resonances,  thus  changing  the  optical  force.  Fig  2.8(a)  shows  the  change  in  force  due  to 
changing  hole  diameter  for  Rl,  e  =  2.28,  slab  thickness  0.6a  and  sphere  radius  of  0.1a. 
The  force  peak  shifts  to  higher  frequencies  as  d  is  increased.  The  force  peak  widens, 
and  the  peak  amplitude  decreases,  corresponding  to  a  decrease  in  the  resonance  quality 
factor.  Fig  2.8(b)  shows  the  change  in  force  with  slab  thickness  for  Rl,  e  =  2.28,  hole 
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Figure  2.7:  (a)  Vertical  force  as  a  function  of  e  for  a  sphere  with  radius  0.15a,  (b)  vertical  force 
as  a  function  of  sphere  radius  for  e=3,  (c)  Vertical  force  for  R1  as  a  function  of  e  for  a  sphere 
radius  of  0.35a. 


diameter  of  0.4a,  and  sphere  radius  of  0.1a. 


2.6  Experimental  considerations 

The  photonie-erystal  slabs  studied  here  are  fabrieated  in  a  silieon-on-insulator  wafer 
using  standard  teehniques  (33).  The  lattiee  eonstant  of  the  photonie  crystal  can  be 
scaled  to  place  a  given  guided  resonance  at  a  particular  wavelength  of  interest  (1).  For 
a  =  760nm,  for  example,  R1  occurs  at  A  =  1609  nm,  while  R2  occurs  at  1532  nm. 
For  a=  375  nm,  R1  is  at  794  nm  and  R2  is  at  748  nm.  An  upper  bound  on  the  power 
required  for  trapping  was  obtained  by  setting  the  height  of  the  potential  depth  |  Af/|  > 
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Figure  2.8:  (a)  Vertical  force  as  a  function  of  e  for  a  sphere  with  radius  0.15o,  (b)  vertical  force 
as  a  function  of  sphere  radius  for  e=3,  (c)  Vertical  force  for  R1  as  a  function  of  e  for  a  sphere 
radius  of  035a. 

10  KbT.  The  power  density  for  a  =  160  nm  was  1  mW/unit  eell  for  Rl,  whieh  has  a 
modest  Q  of  ~260.  The  required  power  deereases  with  inereasing  Q.  Moreover,  the 
power  ealeulated  in  this  manner  is  likely  to  be  overestimated.  Our  results  indieate  that 
the  partieles  will  be  attraeted  to  the  surfaee  of  the  slab.  As  seen  in  Fig  2.3,  the  intensity 
inereases  eloser  to  the  slab  surfaee,  produeing  a  higher  foree  than  the  values  shown  in 
Fig  2.4,  whieh  were  ealeulated  for  a  height  of  0.25a  (for  a  partiele  touehing  the  slab 
surfaee,  with  a  height  of  0.1a). 


In  eonelusion,  we  prediet  that  optieal  forees  above  photonie  erystal  slabs  will  lead 
to  a  variety  of  eomplex,  stably  trapped  erystal  patterns.  Changing  the  wavelength  or 
polarization  of  the  ineident  light  may  be  used  to  reeonfigure  the  patterns.  We  expeet  the 
approaeh  of  light-assisted  self-assembly,  deseribed  here,  to  be  of  broad  utility  for  fab- 
rieation  of  eomplex  photonie  materials,  sensors,  and  filters.  It  is  intriguing  to  eonsider 
whether  metamaterials  based  on  metal  nanoelusters  (34),  for  example,  eould  be  assem¬ 
bled  in  the  manner  we  deseribe,  using  the  optieal  response  of  the  nanoelusters  to  tune  or 
adjust  assembly. 
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Chapter  3 

Experimental  demonstration  of  Light 
Assisted  Template  Self-assembly 


In  this  chapter,  we  demonstrate  for  the  first  time  the  trapping  of  multiple  partieles  using 
an  extended  mode  in  a  two-dimensional  photonie  erystal  slab.  Our  system  improves 
upon  previous  optieal  trapping  work  in  three  main  ways.  First,  we  have  demonstrated 
a  eompaet,  multipartiele  trapping  system  in  the  present  work;  over  180  partieles  are 
trapped  in  an  area  less  than  150  /rm^.  The  number  of  partieles  eould  be  increased  in  a 
straightforward  way  by  proper  design  of  the  strueture  or  by  inereasing  the  size  of  the 
ineident  beam.  Seeond,  since  the  large  array  of  particles  is  trapped  by  exeiting  a  single 
mode  of  the  strueture  without  the  use  of  a  spatial  light  modulator  or  any  movement  of 
the  beam  or  sample,  it  ean  be  easily  adapted  for  a  variety  of  integrated,  lab-on-a-ehip 
applieations.  Finally,  by  tuning  the  polarization  of  the  ineident  light,  we  ean  tune  the 
preferential  direetion  of  the  in-plane  trapping  stiffness  of  our  structure. 

The  proeess  of  light-assisted  templated  self-assembly  (LATS)  is  shown  sehemat- 
ically  in  Fig  3.1.  Light  is  incident  from  below  on  a  photonic-erystal  slab.  The 
photonie-erystal  slab  eonsists  of  a  silieon  deviee  layer  patterned  with  a  periodie  array 
of  air  holes.  The  slab  supports  guided-resonanee  modes,  eleetromagnetie  modes  for 
whieh  the  light  intensity  near  the  slab  is  strongly  enhaneed.  In  previous  work  (35;  36), 
we  have  predieted  theoretieally  that  when  the  incident  laser  is  tuned  to  the  wavelength 
of  a  guided-resonanee  mode,  nanoparticles  will  be  attraeted  toward  the  slab.  The 
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Figure  3.1:  Schematic  picture  of  templated,  light-assisted  self  assembly.  Incident  light  from 
below  excites  a  guided-resonance  mode  of  a  photonic-crystal  slab,  giving  rise  to  optical  forces 
on  nanoparticles  in  a  solution.  Under  the  influence  of  the  forces,  the  nanoparticles  self  assemble 
into  regular  patterns. 

attractive,  optical  force  arises  from  a  strong  electric  field  gradient  just  above  the  slab 
surface.  In  addition,  the  nanoparticles  will  experience  lateral  optical  forces  due  to 
the  electromagnetic  field  structure  of  the  guided  mode.  We  expect  the  optical  forces 
to  result  in  the  assembly  of  particles  into  regular  2D  arrays.  Importantly,  the  array 
patterns  formed  can  be  different  from  the  triangular,  close-packed  structure  formed  by 
traditional,  colloidal,  self  assembly.  Figure  1  depicts  a  square  array  of  particles,  as  one 
example.  By  exciting  different  resonant  modes  of  the  slab,  or  by  designing  different 
slab  templates,  we  expect  that  a  variety  of  patterns  can  be  produced.  Moreover,  these 
patterns  can  reversibly  be  assembled  or  disassembled,  suggesting  the  possibility  of 
’’programmable  optical  matter”. 


3.1  Design 

We  designed,  fabricated,  and  characterized  a  photonic-crystal  slab  for  use  in  the  LATS 
process.  An  electron  micrograph  is  shown  in  Fig  3.2(a).  The  device  was  fabricated 
in  silicon  using  electron-beam  lithography  and  reactive  ion  etching  (see  Methods 
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Figure  3.2:  (a)  SEM  image  of  photonic-crystal  slab.  The  scale  bar  in  the  inset  is  1  /rm.  (b)  3D 
FDTD  simulation  of  the  magnetic  field  {Hz  )  for  a  normally-incident,  x-polarized  plane  wave. 
Circles  represent  the  positions  of  holes;  four  unit  cells  are  shown,  (c)  Measured  transmission 
spectrum  (log  scale). 


in  Sec.  3.4).  The  dimension  and  spacing  of  the  holes  were  designed  to  support 
doubly-degenerate  guided-resonance  modes  near  1.55  /im.  The  magnetic  field  profile 
resulting  from  an  x-polarized,  incident  plane  wave  is  shown  in  Fig  3.2(b).  Fields  were 
calculated  using  the  three-dimensional  finite-difference  time-domain  method  (FDTD). 
The  field-profile  for  a  y-polarized  incident  wave  is  rotated  by  90  degrees.  Fig  3.2(c) 
shows  the  measured  transmission  spectrum  of  the  device.  The  guided-resonance  mode 
appears  as  a  dip  in  the  spectrum.  The  quality  factor,  Q,  was  determined  to  be  ~  170  by 
fitting  to  a  Fano-resonance  shape. 


3.2  Results 

We  carried  out  the  assembly  process  in  a  microfluidic  chamber  filled  with  520-nm 
diameter  polystyrene  particles,  using  a  laser  power  of  64  mW.  Fig  3.3  shows  snapshots 
of  the  DATS  process.  The  photonic-crystal  lattice  is  visible  in  the  background  of  each 
frame.  When  the  laser  beam  is  turned  on,  nanoparticles  are  attracted  towards  the  slab, 
and  begin  to  occupy  sites  within  the  square  lattice  (Fig  3.3(a)).  As  time  progresses. 
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a)  b)  c)  d) 


Figure  3.3:  Light-assisted,  templated  self  assembly  of  520  nm  diameter  partieles  above  a 
photonie-erystal  slab.  The  square  lattiee  of  the  slab  is  visible  in  the  baekground,  oriented  at 
45  with  respeet  to  the  eamera.  (a-e)  Sequential  snapshots  taken  with  the  light  beam  on.  (d) 
Snapshot  taken  after  the  beam  is  turned  off. 

additional  particles  diffuse  into  the  region  where  the  beam  intensity  is  high,  and  begin 
to  form  a  cluster  (Fig  3.3(b)).  Eventually,  a  regular  array  of  partieles  is  formed  (Fig 
3.3(e)).  The  square  symmetry  of  the  assembled  partieles  is  evident  from  the  pieture. 
When  the  laser  beam  is  turned  off,  the  partieles  immediately  begin  to  disperse  and 
diffuse  away  from  the  slab  (Fig  3.3(d)).  The  frames  in  Fig  3.3  were  reeorded  with 
a  dilute  partiele  solution  for  elarity  of  imaging,  and  represent  an  elapsed  time  of 
approximately  one  hour.  Faster  eluster  formation  oeeurs  within  solutions  of  higher 
partiele  eoneentration. 

Each  site  of  the  square  lattiee  may  be  viewed  as  an  optieal  trap.  We  used  partiele- 
traeking  software  to  analyze  partiele  motion  for  fully-assembled  elusters  (see  See.  3.4). 
Fig  3.4  shows  the  reeorded  partiele  positions  extraeted  from  a  20-seeond  video.  The 
ineident  light  was  polarized  along  the  x-direetion  of  the  lattiee.  The  figure  shows  that 
the  partieles  tend  to  stay  above  the  holes  in  the  photonie  erystal,  with  some  variation  in 
position  over  time.  Eaeh  blue  ellipse  represents  a  fit  to  the  data  in  a  single  unit  eell.  It 
ean  be  seen  that  the  variation  in  partiele  position  inereases  at  the  edge  of  the  trapping 
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a)  b)  c) 


Figure  3.4:  Trap  stiffness  for  incident,  x-polarized  light,  (a)  Particle  positions  (red  dots) 
extracted  from  a  20-second  video.  Blue  ellipses  represent  two  standard  deviations  in  position. 

(b)  Histogram  of  stiffness  values  in  the  direction  parallel  to  the  polarization  of  the  incident  light. 

(c)  Stiffness  in  the  direction  perpendicular  to  the  incident  polarization. 
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0.29 

Table  3. 1 :  Polarization  dependence  of  trap  stiffness. 


region,  due  to  the  reduetion  in  power  away  from  the  eenter  of  the  beam. 


The  stiffness  of  eaeh  trap  can  be  determined  from  the  variance  in  particle  position 
(37).  Figs  3.4(b)  and  3.4(c)  show  histograms  of  the  in-plane  stiffness  values  extracted 
from  the  videos.  The  stiffness  of  each  trap  was  normalized  to  the  local  intensity  in  that 
unit  cell  (see  Sec.  3.4).  We  observe  that  the  power-normalized  stiffness  over  the  array 
of  traps  is  normally  distributed,  both  for  the  parallel  and  perpendicular  stiffness.  The 
mean  parallel  stiffness  is  lower  than  the  perpendicular  stiffness,  as  shown  in  Table  3.1 
(O''  angle). 
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We  observe  that  the  trap  stiffness  can  be  tuned  by  rotating  the  direction  of  incident 
light.  When  the  incident  light  is  polarized  at  45°  with  respect  to  the  lattice  directions,  the 
stiffness  values  are  approximately  the  same  in  the  x-  and  y-directions  (Table  3.1).  This 
is  to  be  expected,  since  the  incident  light  excites  both  of  the  doubly-degenerate  modes 
with  equal  strength.  At  90°,  the  stiffness  in  the  perpendicular  direction  (y)  is  again 
lower  than  in  the  parallel  (x)  direction.  The  ability  to  tune  the  stiffness  with  incident 
light  indicates  the  strong  optical  nature  of  our  traps.  The  mean  values  of  stiffness 
are  comparable  to  those  reported  elsewhere  in  the  literature  for  single  particle  traps  (38). 

Using  the  Stokes  drag  method  (see  Sec.  3.4),  we  experimentally  estimated  the 
maximum  force  exerted  on  the  particles  by  the  traps  to  be  0.3  pN. 


3.3  Discussion 

To  understand  how  the  optical  forces  result  in  the  observed  nanoparticle  patterns,  we 
calculated  the  force  numerically  (see  Sec.  3.4).  Fig  3.5(a)  shows  the  force  on  a  520  nm 
diameter  particle  whose  bottom  edge  is  25  nm  above  the  surface  of  the  photonic  crystal 
slab  for  x-polarized  light.  The  background  color  represents  the  vertical  force,  where 
a  negative  force  indicates  attraction  toward  the  slab.  There  are  two  regions  above  and 
below  the  hole  where  the  force  is  slightly  repulsive,  but  at  any  other  position  within  the 
unit  cell  the  particle  results  in  attraction  towards  the  slab.  The  arrows  represent  in-plane 
forces. 

To  determine  the  equilibrium  position  of  the  trapped  particle,  we  calculate  the 
optical  potential.  Given  the  relative  size  of  the  particles  (260  nm  radius)  and  holes  (150 
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a)  b)  c) 


Figure  3.5:  (a)  Force  map  in  one  unit  cell.  The  black  circle  represents  the  position  of  a  hole, 
and  the  color  bar  represents  the  vertical  optical  force  in  dimensionless  units  of  Fc/P,  where  c 
is  the  speed  of  light,  and  P  is  the  incident  optical  power.  A  negative  force  indicates  attraction 
toward  the  slab.  The  arrows  represent  the  magnitude  of  the  lateral  optical  force.  The  length  of 
the  white  arrow  at  the  bottom  center  of  the  figure  corresponds  to  the  maximum  in-plane  value 
of  1.6.  (b)  The  potential  through  the  center  of  the  hole  along  the  x-direction.  The  green  dashed 
lines  indicate  the  position  of  the  hole,  (c)  The  potential  through  the  center  of  the  hole  along  the 
y-direction. 


nm  radius)  in  our  experiment,  the  partiele  can  be  drawn  into  the  hole  by  the  attractive 
vertical  force.  We  calculate  the  optical  potential  as  a  function  of  x-y  position.  For 
each  x-y  position,  the  vertical  height  of  the  particle  is  as  small  as  possible,  given  the 
geometrical  constraints  (see  inset  in  Fig  3.5(b)).  The  result  is  shown  by  the  blue  line 
(total  potential)  in  Figs  3.5(b)  and  3.5(c).  It  can  be  seen  that  the  stable  equilibrium 
position  is  at  the  center  of  the  hole  (x  =  0,  y  =  0),  in  agreement  with  experimental 
observations. 


The  ability  of  the  particle  to  sink  into  the  hole  is  a  key  factor  in  determining  the 
equilibrium  positions.  For  comparison,  the  red,  dashed  line  (in-plane  potential)  in  Fig 
3.5(b)  shows  the  optical  potential  calculated  at  a  constant  z  height  (bottom  edge  of 
particle  25  nm  above  the  slab  surface).  Two  local  minima  are  observed  at  the  edges  of 
the  hole,  which  are  indicated  by  green,  dashed  lines.  From  inspection  of  Fig  3.5(a),  it 
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can  be  seen  that  these  points  correspond  to  locations  where  the  in-plane  forces  are  zero. 
However,  these  two  minima  are  not  stable  equilibrium  positions  of  the  total  potential, 
(blue  line.  Fig  3.5(b-c)). 

From  our  experiments,  we  determined  that  the  threshold  intensity  for  trapping  was 
134  /iW  per  unit  cell  (see  Sec.  3.4).  For  this  intensity,  the  calculated  potential  depth  is 
4.5  keT. 


3.4  Methods 

3.4.1  Fabrication 

The  photonic-crystal  device  consists  of  a  square  lattice  of  holes  etched  into  a  SOI  wafer. 
The  lattice  constant  a  is  860  nm,  and  the  hole  radius  is  0.174a  (150  nm).  The  thickness 
of  the  silicon  device  layer  is  250  nm.  The  buried  oxide  layer  is  3  /rm  thick,  and  the 
silicon  substrate  is  600  /im  thick.  To  fabricate  the  device,  a  layer  of  PMMA  4%  in 
anisole  was  spin  coated  onto  an  SOI  wafer.  Electron  beam  lithography  at  30  kV  was 
used  to  write  the  photonic-crystal  pattern  into  the  resist.  A  modified  Bosch  process  was 
used  to  transfer  the  pattern  into  the  Si  device  layer  using  a  mixture  of  SF6  and  C4F8 
gases.  Plasma-enhanced  chemical  vapor  deposition  was  used  to  deposit  a  193  nm  layer 
of  SiNx  onto  the  polished,  back  surface  of  the  sample  to  reduce  unwanted  reflections. 

The  photonic-crystal  sample  is  mounted  on  a  circular  glass  slide  and  inserted  in  a 
rotary  stage.  The  sample  is  covered  with  a  solution  of  520  nm  diameter  polystyrene 
particles  (700  /xL  of  Thermo  Scientific  Fluoro-Max  R520  particles  in  suspension, 
diluted  in  60  mL  of  deionized  water).  600  /iL  of  1%  Triton-X  was  added  to  the 
solution  to  minimize  particle  stiction  to  the  sample  surface.  On  a  second  glass  slide,  an 
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open-topped  microfluidic  chamber  (4  mm  x  4  mm)  was  fabricated  in  a  5  jim  layer  of 
PDMS  using  photolithography.  The  PDMS  chamber  was  pressed  on  to  the  sample  and 
sealed  inside  the  rotary  stage. 


3.4.2  Optical  setup 

A  Santee  tunable  laser  with  a  tuning  range  from  1500  nm  to  1620  nm  was  connected 
to  an  erbium-doped  fiber  amplifier  and  a  tunable  bandpass  filter  with  a  width  of  1  nm. 
An  adjustable,  neutral-density  filter  and  polarization-control  optics  were  used  to  control 
the  power  and  polarization  of  the  beam.  A  fiber-to-free-space  collimator  directed  the 
beam  to  the  entrance  aperture  of  a  20X  objective  (NA  =  0.26),  which  focused  the  beam 
onto  the  top  surface  of  the  sample  from  the  back  side.  A  second  20X  objective  was 
used  in  conjunction  with  a  beam  splitter  to  collect  the  light  from  the  top  surface  and  to 
image  the  particle  motion  on  a  CMOS  camera.  Prior  to  each  self-assembly  experiment, 
the  transmission  was  measured  in  cross-polarization  mode.  Polarizers  before  and  after 
the  sample  were  oriented  at  90  degrees  from  one  another,  and  the  wavelength  of  the 
transmission  peak  was  identified.  This  laser  was  then  tuned  to  the  peak  wavelength  to 
carry  out  light-assisted  self  assembly. 


3.4.3  Stiffness  analysis 

After  the  assembly  of  a  cluster,  we  recorded  videos  with  a  fixed  exposure  time  of  33  ms, 
the  fastest  value  available  in  our  experimental  set-up.  Typical  videos  were  600  or  more 
frames  in  length.  Particle  motion  was  analyzed  using  MATLAB  algorithms  written 
by  Blair  and  Dufresne  (http://physics.georgetown.edu/matlab/),  which  are  modified 
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versions  of  IDL  routines  written  by  Crocker  and  Grier  (39). 

The  blue  ellipses  in  Fig  LATS2-4(a)  are  obtained  by  fitting  the  data  in  each  unit  cell. 
The  direction  and  relative  lengths  of  the  major  and  minor  axis  are  determined  from  the 
eigenvectors  and  eigenvalues  of  the  scaled  covariance  matrix.  The  ellipse  is  drawn  to 
represent  a  95%  confidence  interval:  positions  outside  the  ellipse  will  be  observed  only 
5%  of  the  time  by  chance  if  the  underlying  distribution  is  Gaussian. 

The  measured  variances  were  corrected  for  motion  blur  due  to  the  finite  integration 
time  of  the  camera  (40).  We  observed  that  the  variance  as  a  function  of  position  within 
the  cluster  had  a  2D  Gaussian  distribution.  This  is  to  be  expected  due  to  the  spatial 
variation  in  intensity  in  the  incident  beam.  We  confirmed  via  3D-FDTD  simulations 
of  a  finite-size  structure  with  Gaussian  beam  excitation  that  the  field  intensity  above 
the  slab  has  a  Gaussian  envelope.  The  stiffness  values  shown  in  Fig  LATS2-4  are 
normalized  to  the  local  power  at  each  trapping  site,  as  determined  by  fitting  the 
experimentally-measured  variances  to  a  2D  spatial  Gaussian. 

We  directly  measured  the  diffusion  coefficient  by  first  assembling  a  cluster,  then 
blocking  the  beam  and  performing  a  linear  regression  fit  to  the  subsequent  diffusion  as 
a  function  of  time  (39).  A  value  of  0.56  fim? / s  was  obtained. 


3.4.4  Force  measurement 

After  assembling  a  cluster,  we  increased  the  flow  speed  until  all  particles  were  released. 
We  observed  the  release  of  the  last,  trapped  particle  at  a  fluid  velocity  of  30  /rm/s  for 
a  laser  power  of  64  mW.  In  the  low  Reynolds  number  regime,  the  particle  velocity  v  is 
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related  to  the  external  force  by  the  equation  F  =  QTirjrv,  where  rj  is  the  viscosity  of  the 
medium,  and  r  is  the  radius  of  the  particle. 


3.4.5  Force  and  potential  calculation 

For  each  particle  position  within  the  unit  cell,  we  performed  an  FDTD  simulation  to 
calculate  the  electromagnetic  fields  and  obtained  the  optical  force  from  an  integral 
of  the  Maxwell  Stress  Tensor  over  a  box  just  large  enough  to  include  the  particle. 
The  z  dependence  of  the  force  inside  the  hole  was  obtained  from  an  exponential  fit 
to  the  field  intensity  above  the  slab.  The  power  threshold  for  stable  trapping  was 
determined  experimentally  from  the  size  of  the  stable  cluster  (16  unit  cells,  or  13.8  /rm, 
in  diameter).  The  Gaussian  distribution  of  energy  in  the  mode,  obtained  from  a  fit  of 
the  variance  measurements,  had  a  FWHM  of  17  /rm.  Given  the  input  power  of  64  mW 
we  find  a  peak  intensity  of  240  iiW  per  unit  cell,  and  a  trapping  threshold  intensity  at 
the  edge  of  the  cluster  of  134  jjW  per  unit  cell. 


3.5  Conclusions 

In  summary,  we  have  demonstrated  the  technique  that  we  term  light-assisted,  templated 
self-assembly.  A  square  array  of  ~  180  polystyrene  particles  was  assembled  above  a 
square-lattice  photonic  crystal.  The  assembly  of  the  array  resulted  from  the  near  field 
of  an  extended  mode  in  a  photonic  crystal  slab.  The  optical  force  was  enhanced  at  the 
resonant  frequency  of  the  mode.  We  characterized  the  trapping  stiffness  of  our  system, 
and  we  showed  that  the  calculation  of  optical  forces  successfully  predicts  the  observed 
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equilibrium  trapping  positions. 

Our  system  can  be  re-designed  to  assemble  larger  clusters.  One  approach  would 
involve  using  a  mode  with  a  higher  quality  factor.  The  optical  force  scales  with  field 
intensity,  and  higher-Q  modes  should  have  higher  local  field  intensity  in  the  near 
field  for  a  given  input  power.  Furthermore,  photons  in  higher-Q  modes  spread  out 
further  within  the  slab  due  to  their  longer  lifetimes,  thereby  increasing  the  trapping 
area.  Another  approach  is  to  re-design  the  photonic-crystal  slab  to  increase  the  field 
concentration,  as  shown  previously  (36).  Expanding  or  re-shaping  the  incident  beam 
may  also  permit  the  formation  of  larger  clusters. 

One  major  advantage  of  our  approach  to  self-assembly  is  the  ability  to  assemble 
complex  structures  with  symmetries  that  are  not  usually  available  due  to  energy 
concerns;  here  we  have  demonstrated  just  one.  Different  photonic-crystal  templates 
will  lead  to  different  particle  patterns  and  more  complex  unit  cells.  For  a  given 
template  the  particle  patterns  can  be  reconfigured  by  changing  the  wavelength  or 
polarization  of  the  incident  beam.  The  nature  of  the  particle  itself  can  also  be  changed; 
interesting  applications  may  result  from  the  use  of  metal  nanoparticles,  quantum  dots 
or  combinations  of  different  species  (35). 

Our  system  also  lends  itself  to  compact  integration  on-chip.  The  only  external 
component  in  our  experiment  is  the  laser  beam,  and  by  fabricating  the  device  in  an 
active  material,  the  source  itself  could  also  be  integrated.  Using  electrical  excitation 
of  such  a  photonic-crystal  laser  would  further  miniaturize  our  structure  in  a  way  that 
would  enable  a  variety  of  interesting  on-chip  applications. 
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We  also  anticipate  that  our  technique  will  find  applications  in  the  fabrication  of 
metamaterials  and  other  photonic  devices.  After  assembling  a  2D  array  of  particles, 
polymerization  could  be  used  to  transfer  the  array  to  another  substrate.  The  technique 
could  thus  be  used  to  create  3D  arrays  layer-by-layer.  With  suitable  design  of  the 
structure  it  should  also  be  possible  to  assemble  3D  arrays  in  situ. 

Finally,  a  variety  of  dynamic,  real-time  applications  could  be  enabled  by  DATS.  For 
instance,  LATS  could  be  used  to  form  reconfigurable  photonic  materials,  such  as  a  tun¬ 
able  photonic  filter.  The  pattern  of  particles  assembled  can  be  altered  in  a  reversible 
way  simply  by  changing  the  input  polarization  or  wavelength.  By  combining  the  opti¬ 
cal  potential  landscape  with  microfluidic  flow  channels,  it  should  also  be  possible  to 
sort  particles  into  different  channels  based  upon  size  or  refractive  index  (18;  19).  The 
manipulation  of  biological  objects  could  be  simplified  in  our  system  by  allowing  for 
batch  processing  of  tasks  which  typically  use  a  single  optical  trap  to  examine  a  sin¬ 
gle  specimen.  Furthermore,  by  trapping  biological  objects  with  controlled  spacing,  our 
system  allows  for  control  over  interactions. 
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Chapter  4 

Optical  trapping  of  Metal-dielectric 
cluster  near  a  micro-cavity 


Optical  trapping  and  manipulation  of  particles  has  been  demonstrated  in  a  large  number 
of  struetures  sinee  the  diseovery  of  optieal  tweezers.  Resonant  struetures  have  been 
of  great  interest  for  trapping,  sorting  and  eooling  of  partieles  with  different  optieal 
properties  sueh  as  atoms  (41),  quantum  dots  (42),  DNA  (43),  resonant  partieles  (44;  45), 
and  dieleetrie  partieles  (46).  Different  meehanism  are  used  to  resonantly  enhanee  the 
eleetromagnetie  field  near  a  structure,  in  partieular,  waveguides  (47),  ring  resonators 
(48),  nanoantennas  (49;  50),  plasmonie  struetures,  photonie  erystals  and  photonie 
erystals  eavities  (51;  52)  are  the  most  popular  ehoiees. 

However,  in  the  dipole  approximation,  the  gradient  foree  negleets  the  field  redistri¬ 
bution  inside  the  optieal  eavity  due  to  the  polarization  of  the  trapped  partiele.  In  this 
ehapter,  I  use  a  FDTD  method  to  solve  Maxwell  equations  and  eompute  the  Maxwell 
stress  tensor  for  a  gold  partiele  on  top  of  a  miero-eavity.  Next,  I  show  a  seeondary  trap 
formation  due  to  the  field  redistribution  around  the  trapped  gold  partiele  and  present  an 
optieal  indueed  metal-dieleetrie  eluster  formation. 
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4.1  Trapping  of  gold  nano-particles 

As  mentioned  before,  in  templated  self-assembly  methods  the  geometrical  relation 
between  template  and  particles  mainly  defines  the  final  structure.  Thus,  cluster  for¬ 
mation  of  particle  with  same  geometry  but  different  composition  cannot  be  controlled 
independently.  The  use  of  optical  forces  for  particle  trapping  allows  for  differential 
response  based  on  composition,  since  dielectric  and  metal  particles  experience  different 
optical  forces. 

Strong  optical  forces  can  be  obtained  from  the  field  enhancement  near  micropho¬ 
tonic  structures.  In  particular,  the  microcavity  mode  of  a  photonic  crystal  can  be  used 
to  create  a  nanoparticle  trap  (9;  10).  In  this  section  we  present  a  controllable  formation 
of  a  metal/dielectric  cluster.  We  propose  to  use  optical  forces  on  particles  in  a  solution 
above  a  photonic-crystal  microcavity  and  the  field  redistribution  due  to  the  presence 
of  a  trapped  particle.  First,  The  photonic-crystal  microcavity  forms  an  optical  trap  for 
gold  nano-particles.  Field  enhancement  near  the  trapped  gold  nanoparticle  then  forms 
secondary  trapping  sites  for  a  pair  of  dielectric  nanoparticles.  The  process  allow  optical 
driven  formation  of  a  well-defined  metal/dielectric  cluster. 

In  contrast  to  previous  work  on  optical  binding  where  clustering  comes  from 
interactions  between  dipoles  (27),  in  this  paper  we  are  considering  nondipolar  colloidal 
particles.  Therefore,  a  trapped  gold  nano-particle  on  a  photonic  crystal  microcavity  acts 
as  an  off-resonance  plasmonic  trap  for  dielectric  particles  instead  of  a  dipole-like  parti¬ 
cle.  While  previous  work  has  examined  light  forces  on  metallic  particles  (53;  54;  55), 
between  metallic  nano-particles  (56),  and  forces  due  to  plasmonic  traps  (45;  57),  optical 
forces  in  a  particle-like  off-resonance  plasmonic  trap  were  not  considered. 
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The  structure  we  consider  is  shown  in  Fig  4.1(a).  In  an  experiment,  light  would 
illuminate  the  photonic  crystal  microcavity  from  below.  The  photonic  crystal  is  formed 
of  a  hexagonal  lattice  of  holes  in  Si  (n  =  3.46).  The  lattice  constant  is  400  nm,  the 
hole  radius  is  140  nm,  and  the  thickness  of  the  silicon  slab  is  200  nm.  The  microcavity 
is  created  by  removing  a  hole  and  displacing  its  two  nearest  neighbors  along  the  x 
direction  outward  by  8  nm  (9).  The  photonic  crystal  rests  on  a  Si02  substrate  (n  = 
1.45),  and  the  holes  in  the  lattice  and  the  region  above  the  crystal  are  filled  with  water 
(n=  1.33). 

The  microcavity  mode  was  simulated  using  the  finite-difference  time-domain 
method  (FDTD).  A  y-polarized  plane  wave  illuminates  the  microcavity  from  below. 
We  take  the  computational  cell  to  be  3.6  /rm  x  3.4  /rm  x  2  /rm  in  size.  PML  boundary 
conditions  were  used.  A  microcavity  resonance  value  of  226.9  THz  and  quality  factor 
Q  ~  60  were  obtained.  The  microcavity  strongly  enhances  the  electric  field  near  the 
missing  hole,  and  the  field  intensity  is  80  times  higher  than  without  the  microcavity 
in  the  midplane  of  the  slab.  The  field  intensity  in  (V/m)'^  is  plotted  in  Fig.  4.1(b) 
for  a  source  power  of  l.SAmW/ fj,rrS.  The  intensity  decays  by  a  factor  of  30  at  a 
height  of  130  nm  above  the  slab  Fig  4.1(c).  At  this  height,  we  have  observed  that 
the  strongest  local  intensity  is  achieved  slightly  off  resonance  (f  =  229.1  THz)  Fig  4.1(d). 

A  gold  nano-particle,  far  from  the  plasmon  resonance,  can  be  seen  as  particles 
with  slow  varying  frequency  dependent  imaginary  dielectric  constant,  when  the  size 
of  the  particle  is  bigger  than  the  mean  free  path  of  the  electrons  in  the  metal  (58;  59). 
However,  when  the  particle  size  is  comparable  to  the  electric  field  variation  of  the 
optical  trap,  the  force  cannot  be  approximated  as  the  sum  of  a  gradient,  scattering  and 
absorption  force  due  to  the  strong  field  redistribution.  Therefore,  in  order  to  calculate 
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Figure  4. 1 :  (a)  Schematic  of  defect  cavity  in  a  hexagonal  lattice  of  holes  with  diameter  d  and 
thickness  h,  (b)  Electric  field  distribution  at  the  center  of  the  slab  at  resonance,  (c)  Electric  field 
disfribufion  af  z=130  nm  above  fhe  slab  af  resonance,  (d)  Elecfric  field  disfribufion  af  z=  130  nm 
above  fhe  slab  slighfly  off  resonance. 

the  optical  force  on  spherical  gold  particles  we  used  full-vectorial  calculations  of 
Maxwell’s  equations  and  the  Maxwell  Stress  Tensor  method  (22). 

In  Fig  4.2,  we  plot  the  normalized  optical  force  in  the  vertical  direction  for  gold 
nanoparticles  with  several  different  radii.  The  bottom  edge  of  each  particle  is  50  nm 
above  the  surface  of  the  slab,  and  each  particle  is  centered  on  the  microcavity.  is  the 
z  component  of  the  optical  force,  c  is  the  speed  of  light  in  vacuum,  and  the  flux  is  the 
total  electromagnetic  flux  over  the  computational  cell.  A  negative  force  indicates  attrac¬ 
tion  toward  the  microcavity,  whereas  a  positive  force  indicates  repulsion.  We  observe 
both  strong  size  and  frequency  dependence  of  the  optical  force  near  the  microcavity 
resonance.  For  the  largest  gold  particle  shown  (r  =  160  nm),  the  force  exhibits  a  Fano- 
resonance-like  shape  and  is  repulsive  (positive  force).  For  the  smallest  gold  particle 
shown  (r  =  80nm),  the  force  exhibits  a  single  dip  and  is  attractive  (negative  force).  For 
particles  of  intermediate  sizes,  the  force  transforms  from  a  Fano  shape  to  a  single  dip. 
Thus,  only  gold  particles  with  radii  smaller  than  160  nm  can  be  trapped.  Furthermore, 
the  force  on  an  80nm  gold  particle  is  an  order  of  magnitude  larger  than  that  felt  by  a 
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Figure  4.2:  Vertical  force  on  gold  nanoparticles  of  different  sizes  and  a  dielectric  particle  with  r 
=  80  nm  and  e  =  1.59. 


dielectric  particle  of  the  same  size. 


Qualitatively,  we  can  understand  the  different  shapes  of  the  curves  in  Fig  4.2  using 
insight  from  the  dipole  approximation.  For  trapping  to  occur,  the  attractive  gradient 
force  must  be  larger  than  the  sum  of  the  repulsive  scattering  and  absorption  forces.  For 
small  radii,  the  electric  field  polarizes  the  particle  such  that  the  gradient,  scattering,  and 
absorption  contributions  to  the  optical  force  can  be  written  as  (54): 

Fgrad=Y^\E\^  (4.1) 

^  sea 

C 

"flwSC abs 

^  abs 

C 

From  the  analytical  expressions  in  Ec  (4.1)  we  plot  the  forces  in  Fig  4.3.  It  can  easily 
be  seen  that  the  gradient  force  on  a  gold  particle  is  attractive  and  has  the  form  of  a  single 
dip  near  resonance,  while  the  scattering  and  absorption  forces  are  repulsive  and  have  a 
Fano-like  shape. 

Looking  at  Fig  4.2,  we  observe  that  for  small  particles  (e.g.  r  =  80  nm),  the  shape  of 
the  optical  force  curve  resembles  that  of  a  gradient  force  term.  For  larger  particles  (e.g. 
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Figure  4.3:  Gradient  (a),  and  scattering  and  absorption  forces  (b)  on  a  gold  particle  with  r  = 
80nm  using  the  dipole  approximation. 


r  =  160  nm),  the  shape  of  the  curve  resembles  that  of  scattering  and  absorption  terms. 
We  infer  that  as  the  particle  size  decreases,  the  gradient  force  dominates  the  scattering 
and  absorption  effects,  leading  to  trapping. 


In  Fig  4.4,  we  plot  the  frequency-dependent  vertical  force  on  a  gold  particle  with  r 
=  80  nm  for  several  different  heights,  assuming  the  particle  is  centered  over  the  micro¬ 
cavity.  The  frequency  at  which  the  optical  force  peaks  shifts  as  a  function  of  particle 
height,  due  to  the  interaction  between  the  gold  particle  and  the  microcavity.  This  effect 
is  similar  to  the  self-induced  trapping  effect  previously  observed  for  dielectric  particles 
(9).  For  a  frequency  of  229.1THz,  the  force  is  attractive  for  all  heights  shown. 


We  fix  the  excitation  frequency  at  229.1THz  and  consider  the  spatial  dependence 
of  the  optical  force.  In  Figs  4.4(b)  and  4.4(c),  we  show  the  normalized  horizontal  and 
vertical  forces  (Fc  /  flux)  on  a  gold  particle  whose  bottom  edge  is  located  50nm  above 
the  slab.  The  color  bars  indicate  the  magnitude,  and  the  arrows  in  Fig.  4.4(b)  indicate 
the  direction.  There  is  a  stable  trapping  position  at  the  center  of  the  microcavity.  The 
trap  is  weakest  along  the  y-direction  and  has  a  stiffness  [10]  ky  =  0.2  pNnm~^W~^. 
We  estimate  a  minimum  trapping  power  of  ImW  by  requiring  the  stability  number  to 
be  at  least  one  (60). 
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Figure  4.4:  (a)  Vertical  force  as  a  function  of  frequency  for  several  different  vertical  positions. 
Normalized  force  in  the  (b)  horizontal  and  (c)  vertical  directions  as  a  function  of  position  in  the 
x-y  plane  for  a  gold  particle  with  r  =  80nm. 

4.2  Cluster  formation 

A  gold  nanoparticle  trapped  at  the  eenter  of  the  mieroeavity  will  redistribute  the  elee- 
tromagnetie  field,  forming  new  optical  trapping  locations.  Figs  4.5(a)  and  4.5(b)  show 
the  eleetrie  field  intensity  130nm  above  the  slab  at  the  fixed  frequeney  of  229.1THz  for 
the  same  souree  power  as  in  Fig  4.1.  The  maximum  value  of  intensity  is  7  times  higher 
than  without  a  gold  partiele.  The  high  intensity  regions  ean  attraet  dieleetrie  partieles. 
Figs  4.5(e)  and  4.5(d)  show  the  horizontal  and  vertical  optical  forces  on  a  dieleetrie 
partiele  with  n  =1.59  and  r  =  80  nm  in  the  presenee  of  the  trapped  gold  nanopartiele, 
indieated  by  the  solid,  yellow  eirele.  The  dielectric  particle  is  stably  trapped  in  the 
position  indieated  by  the  dashed  eireles.  In  the  Figures,  the  arrows  next  to  the  graphs 
represent  values  of  Fe  /  flux  =  0.002.  We  estimate  a  minimum  trapping  power  of  7  mW 
with  trapping  stiffnesses  of  =  0.08  pNnm~^W~^ ,  ky  =  0.25  pNnm~^W~^  and  kz= 
0.09  pNnm~^W~^ . 
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Figure  4.5:  Electric  field  intensity  in  the  midplane  (a)  and  on  a  vertical  cut  (b)  of  a  trapped  gold 
particle  with  r  =  80  nm.  Horizontal  (c)  and  vertical  (d)  optical  forces  as  a  function  of  position 
for  a  dielectric  particle  with  r  =  80  nm  in  the  presence  of  a  trapped  gold  particle  (yellow  circle), 
(e)  Horizontal  optical  force  as  a  function  of  position  for  a  dielectric  particle  near  trapped  gold 
(yellow)  and  dielectric  (blue)  particles.  The  arrows  next  to  the  graph  represent  Fc  /  flux  =  0.002. 

By  symmetry,  a  second  possible  trapping  position  exists  on  the  opposite  side  of  the 
gold  particle.  Fig.  4.5(e)  shows  the  optical  force  on  a  dielectric  particle  with  r  =  80  nm 
after  both  a  gold  and  dielectric  particle  are  trapped  in  the  locations  shown  by  the  solid 
yellow  and  blue  circles,  respectively.  The  arrow  next  to  the  graph  represents  a  value 
of  Fc  /  flux  =  0.002.  The  dielectric  particle  is  trapped  in  the  position  indicated  by  the 
dashed  line.  The  results  of  Fig  4.5  thus  indicate  that  dielectric -metal-dielectric  clusters 
will  be  formed,  consisting  of  three  particles  with  composition  and  orientation  as  shown. 


We  note  that  the  specific  force  values  and  the  required  power  for  trapping  calculated 
here  assume  normally-incident  light  distributed  over  an  area  of  approximately  12  /im^. 
In  an  experiment,  the  coupling  efficiency  between  the  incident  light  and  the  cavity  can 
likely  be  increased  by  adjusting  the  spot  size  of  the  normally-incident  beam  and/or 
using  a  side-coupled  waveguide  to  inject  light  into  the  microcavity  mode,  increasing 
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the  optical  forces  and  reducing  the  power  required  for  trapping. 

In  an  experiment,  surface  treatment  of  the  gold  nanoparticles,  dielectric  nanopar¬ 
ticles,  and  photonic  crystal  may  be  required  to  prevent  aggregation,  adhesion,  and/or 
electrostatic  effects.  Previous  experimental  work  in  the  literature  has  developed 
techniques  for  preventing  aggregation  of  gold  (55)  and  dielectric  (61)  nanoparticles,  as 
well  as  adhesion  between  nanoparticles  and  microphotonic  devices  (51).  Electrostatic 
interactions  in  the  system  can  be  minimized  using  a  phosphate  buffer  solution  with  a 
regulated  pH  of  7.0  (47). 

We  estimate  that  the  rms  variation  of  the  gold  particle  due  to  Brownian  motion 
is  approximately  20nm  for  the  trap  stiffness  and  power  given  in  the  paper.  We  have 
calculated  the  field  profile  for  displacements  of  the  gold  particle  on  this  order.  We 
observe  that  the  variation  in  field  intensity  is  less  than  10%.  While  the  symmetry  of 
the  intensity  distribution  is  broken  slightly  due  to  particle  motion,  both  the  qualitative 
shape  and  magnitude  is  largely  preserved.  We  therefore  do  not  expect  Brownian  motion 
to  significantly  affect  the  conditions  for  optical  trapping. 

We  have  also  checked  that  the  force  exerted  on  the  gold  particle  by  the  dielectric 
particle  is  much  smaller  than  the  force  exerted  on  the  gold  particle  by  the  cavity. 
The  Figure  below  shows  the  optical  forces  on  the  gold  (yellow)  and  dielectric  (blue) 
particles  in  the  configuration  at  which  the  forces  are  strongest.  The  arrow  next  to 
the  graph  indicates  Fc/fiux=0.002.  We  calculate  that  the  force  on  the  gold  particle 
will  displace  the  particle  by  less  than  12nm  from  the  center  of  the  photonic-crystal 
microcavity.  As  seen  in  our  answer  to  point  4  above,  displacements  on  this  order  will 
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Figure  4.6:  Force  on  gold  particle  due  to  dielectric  particle, 
not  significantly  change  the  field  profile. 

We  finally,  simulated  the  heating-indueed  fluid  eonvection  in  COMSOL  and  used 
Stokes  law  to  estimate  the  drag  force  exerted  on  particles  by  flow  in  a  10  /im  thiek  ehan- 
nel.  The  heating  source  in  the  simulation  is  ealeulated  using  the  spatially-dependent 
eleetrie  fields  obtained  from  FDTD  simulations.  For  an  ineident  power  of  10  mW,  the 
maximum  temperature  rise  is  ~  7  K,  and  the  maximum  flow  speed  is  ~  0.1  nm/s.  The 
normalized  drag  force  (F^c  /  flux)  on  a  80  nm  radius  gold  partiele  is  estimated  to  be  less 
than  4  x  10“®,  whieh  is  negligible  eompared  to  the  optical  trapping  force. 

In  summary,  we  have  proposed  a  method  for  using  optieal  forees  to  induee  the  for¬ 
mation  of  metallo-dieleetrie  elusters  with  eontrolled  morphology.  The  method  relies  on 
field-enhaneement  effeets  near  gold  partieles  trapped  above  a  photonie-erystal  miero- 
eavity.  If  the  self-assembly  proeess  is  performed  in  a  photopolymerizable  solution,  a 
seeondary  laser  spot  could  be  used  to  loeally  solidify  the  region  around  the  eluster. 
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Moreover,  integration  with  a  microfluidic  channel  could  allow  for  controlled,  sequen¬ 
tial  delivery  of  individual  gold  and  metal  particles,  along  with  controlled  collection  of 
metallo-dielectric  clusters. 
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Chapter  5 

GUV’s  Bending  modulus  measurement 
using  Dual-beam  optical  trap 


Dual  beam  optical  trap  (DBOT)  has  seen  a  growing  interest  in  the  field  of  optical  manip¬ 
ulation  of  biological  cells  due  to  its  ease  of  implementation  and  non-invasive  nature. 
Cells  (70)  and  giant  unilamillar  vesicles  (GUVs)(71),  have  been  successfully  trapped 
and  stretched  in  a  DBOT  platform,  which  provide  an  excellent  tool  to  investigate  the 
membrane  elastic  properties.  Thus,  the  DBOT  provides  a  method  for  non-invasive 
application  of  time-dependent  forces  in  a  device  suitable  to  rapid,  high-throughput 
measurements. 

GUVs  typically  contain  the  same  solution  on  the  interior  and  exterior  of  the  spheri¬ 
cal,  lipid  bilayer,  thus  the  refractive  index  difference  required  for  trapping  is  not  present. 
Here,  we  used  GUVs  with  different,  osmotically -balanced  solutions  on  the  interior  and 
exterior  of  the  bilayer.  We  construct  a  DBOT  based  on  optical  fibers  integrated  with  a 
capillary  flow  channel  and  demonstrate  trapping  and  stretching  of  a  GUV.  We  charac¬ 
terize  the  response  of  the  GUV  to  a  step  increase  in  stress.  We  further  conduct  exper¬ 
iments  in  which  the  area  strain  of  the  GUV  is  measured  as  a  function  of  applied  stress 
and  develop  a  method  to  extract  the  bending  modulus  of  the  membrane  from  the  data. 
This  method  incorporates  three-dimensional  ray-tracing  methods  to  calculate  the  applied 
stress  in  the  DBOT.  The  value  we  obtain  for  the  bending  modulus  agrees  well  with  liter¬ 
ature  values.  Unlike  previous  results  in  the  literature  (62;  63),  we  use  a  laser-wavelength 
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Beam  trajectory 


Figure  5.1:  Dual-beam  optical  trap 

selected  to  minimize  heating  effects  (64;  65)  avoid  the  use  of  polyethylene  glycol,  which 
is  known  to  destabilize  membranes  (66),  and  model  the  mechanical  properties  of  the 
GUV  as  a  lipid  bilayer  according  to  the  accepted  approach  developed  by  Helfrich  (67). 
Our  results  demonstrate  the  potential  of  the  DBOT  for  rapid,  flow-through  measure¬ 
ments  of  membrane  response  to  changing  physiochemical  environments,  opening  a  path 
for  a  wide  range  of  biological  experiments. 

5.1  Stress  profile 

Counter-propagating  beams  emitted  from  two  optical  fibers  form  an  optical  trap  due 
to  their  Gaussian  intensity  profile  and  radiation  pressure  (Fig  5.1)(72).  Though  the 
net  force  acting  on  the  object  in  this  configuration  is  zero,  it  experiences  a  non-zero 
surface  stress  (73).  This  surface  stress  deforms  the  surface  by  an  amount  depending 
upon  intrinsic  elasticity.  Thus,  elastic  objects  like  cells  and  GUVs  assume  the  shape  of 
a  prolate  spheroid  (67)  when  trapped  with  an  eccentricity  dependent  upon  the  trapping 
power  and  elasticity.  In  order  to  quantify  the  shape  deformation  of  GUVs  and  cells,  it 
is  of  utmost  importance  to  calculate  the  surface  stress  profile  on  spheroidal  objects  with 
a  refractive  index  comparable  to  the  cells  and  GUVs.  An  exact  calculation  of  optical 
stress  profile  on  the  surface  of  spheroidal  objects  is  an  absolute  requirement  to  quantify 
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the  resultant  shape  deformation  of  the  object  based  on  its  intrinsic  elasticity. 


5.1.1  Surface  stress  using  Ray  Optics 

We  consider  a  spheroidal  dielectric  object  illuminated  by  a  single  Gaussian  beam.  The 
spheroidal  surface  is  resolved  into  grid  points  with  corresponding  azimuthal  angle  6 
varying  from  0  to  7r/2  and  polar  angle  0  from  0  to  27r.  In  the  case  of  a  closed  object, 
a  transmitted  beam  experiences  multiple  internal  reflections.  If  we  consider  m  internal 
reflections,  the  total  optical  surface  stress  is  (29) 
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The  incident  vector  takes  in  to  account  the  divergence  of  the  beam,  meanwhile  the 
subsequent  reflected  and  transmitted  vectors  can  be  computed  using  Snell’s  law.  The 
stress  due  to  the  first  reflection  from  the  incident  beam  is  computed  and  subsequent 
surface  stress  arising  from  multiple  reflection  of  the  transmitted  beam  are  also  recorded. 
Given  the  symmetries  of  the  system  and  provided  non-polarized  light,  only  the  stress 
calculations  in  the  equatorial  plain  are  needed,  any  other  cut  of  the  spheroid  would  be  a 
rotation  of  the  equatorial  plane. 


Fig  5.2,  shows  the  optical  stress  from  the  incident  beam  (ai)  and  the  first  four 
internal  reflection  ((T2  —  for  a  sheroid  of  volume  equal  to  523.6  prc?  and  e  =  0.4. 
The  angular  range  [0  —  7r/2]  corresponds  to  the  illuminated  part  of  the  object  while 
[7r/2  —  tt]  is  the  shadowed  region.  For  ai  there  is  a  maximum  stress  in  the  shadowed 
region  since  a  larger  change  in  momentum  is  observed  when  a  ray  goes  from  a  higher 
to  lower  refractive  index.  Next,  the  surface  stress  from  each  bounce  is  added  vectorialy 
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Figure  5.2:  (a)  Optical  stress  for  multiple  reflections  on  a  spheroid  with  volume  equal  to  a  5  /rm 
sphere  and  e  =  0.4  (ni  =  1.335,  n2  =  1.38,  w  =  \  \  pm).  RI  of  Red  Blood  Cells  is  close  to  the 
value  chosen  in  this  case,  (b)  Total  optical  surface  stress  in  the  equatorial  plane  for  single  beam, 
and  (c)  total  optical  surface  stress  for  dual  beam. 

to  obtain  the  total  optical  surface  stress  profile  in  the  equatorial  plane  from  a  single 
beam  (Fig  5.2(b)).  As  we  can  see,  each  reflection  introduces  a  spike  in  the  total  surface 
stress  near  the  region  of  illumination.  These  spikes  can  be  understood  as  an  interference 
phenomenon  inside  the  spheroid  for  each  single  beam. 

In  an  experiment,  a  low  coherence  source  is  useful  for  dual  beam  optical  trapping 
since  it  eliminates  interference  effects.  The  two  trapping  beams  are  incoherent  with 
respect  to  each  other,  since  the  path  length  difference  between  the  two  fibers  are  much 
longer  than  the  coherence  length  of  the  light.  Thus,  the  stress  contribution  from  the 
dual  beam  trap  is  obtained  via  vectorial  addition  of  the  stress  due  to  the  two  opposing 
Gaussian  beams  (Fig  5.2(c))).  As  expected,  the  optical  stress  is  symmetric  and  the 
spikes  represent  interferences  from  each  beam,  but  not  between  the  beams.  In  order 
to  corroborate  this  result  we  performed  the  same  calculation  using  the  Maxwell  Stress 
tensor  method. 
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5.1.2  Surface  stress  using  the  Maxwell  Stress  tensor 

The  Maxwell  Stress  tensor  method  (MST)  express  the  optical  force  acting  on  an  object 
via  the  change  in  electromagnetic  momentum  across  the  surface  of  a  close  object.  If  we 
want  to  consider  the  optical  stress  in  the  surface,  we  can  consider  two  surfaces,  one  just 
above  the  object  interface  and  one  just  bellow  (28).  The  total  force  is  then 
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Expanding  the  stress  tensor  matrix  in  terms  of  the  associated  electric  and  magnetic  field 
components  and  applying  continuity  of  fields  across  the  surface,  the  time  average  surface 
stress  takes  the  following  form 
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Dn  and  Et  are  the  normal  component  of  electric  flux  density  and  tangential  component 
of  the  electric  field  respectively  and  eo  is  the  permittivity  of  vacuum.  In  agreement  with 
the  RO  model,  the  surface  stress  at  any  point  P  is  directed  along  the  surface  normal  h. 


In  our  simulations,  we  performed  a  complete  full  wave  solution  of  the  Maxwell’s 
equations  using  FDTD  analysis  to  compute  the  electric  fields.  We  simulated  a  circularly 
polarized  light  to  emulate  time-averaged  unpolarized  light.  We  first  simulated  the  single 
beam  case  and  computed  the  corresponding  stress  components.  Using  mirror  symmetry 
we  flipped  the  electric  field  to  simulate  a  counter  propagating  beam  and  calculated  the 
stress  components.  A  scalar  addition  of  these  two  stress  components  (since  both  are 
directed  along  the  surface  normal)  yields  the  total  contribution  from  two  incoherent 
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Figure  5.3:  Circularly  polarized  electric  field  intensity  and  corresponding  stress  profile  in  the 
equatorial  plane  for  (a)  Single  Beam  (b)  Dual  beam. 

counter  propagating  beams. 

Fig  5.3(a)  shows  the  stress  profile  and  the  eorresponding  field  intensity  profile  on 
the  equatorial  plane  of  a  spheroid  objeet  with  refraetive  index  n2  =  1.38  arising  from 
a  single  Gaussian  beam  propagating  in  a  medium  with  refraetive  index  ni  =  1.335 
(initial  beam  waist  wq  =  2.4/rm  and  propagation  distanee  of  110/im  ).  The  dual  beam 
stress  profile  and  the  eleetrie  field  intensity,  obtained  by  ineoherent  sum  of  eleetrie  field 
intensity  arising  from  eaeh  beam  (Fig  5.3(b)). 

5.1.3  Comparison  between  RO  and  MST 

We  earned  out  a  eomparison  in  order  to  verify  the  stress  profile  obtained  using  RO 
and  MST.  We  varied  the  eeeentrieity,  size  of  the  spheroidal  objeet  and  refraetive  index 
eontrast  between  the  surrounding  medium  and  the  objeet.  In  eaeh  of  the  following 
eases,  wavelength  is  set  at  808  nm  and  the  initial  beam  waist  is  wq  =  2A  /rm  and  the 
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Figure  5.4:  Comparison  of  stress  profile  arising  from  single  beam  in  the  equatorial  plane 
obtained  from  ray  optics  and  MST  for  varying  eccentricities  (b)  single  beam  electric  field  inten¬ 
sity  profile  corresponding  to  e  =  0.0  and  e  =  0.9  (c)  Dual  beam. 


medium  is  assumed  to  be  water  which  possesses  a  refractive  index  n  =  1.335.  The 
choice  of  beam  waist  and  wavelength  agree  with  our  experimental  set  up  .  In  case  of 
refractive  index  and  eccentricity  variation  the  trap  center  is  at  distance  of  110/rm  from 
the  beam  waist,  which  makes  the  ratio  of  the  beam  waist,  at  the  trap  site,  and  radius 
of  an  equal  volume  sphere  close  to  1.  In  this  way,  we  establish  similar  illumination 
conditions  for  all  the  cases  studied.  In  first  two  cases,  the  volume  of  the  object  is  kept 
equal  to  that  of  a  sphere  of  radius  5  /rm.  In  case  of  the  size  variation,  the  beam  width  at 
the  trap  center  is  adjusted  by  varying  the  location  of  the  trap  site  to  give  a  ratio  of  beam 
waist  to  object  radius  equal  to  1. 


Fig  5.4  shows  a  comparison  between  RO  and  MST  with  spheroids  of  different 
eccentricities  and  a  refractive  index  contrast  of  0.0105  (similar  to  what  exits  in  case  of 
GUVs  with  sucrose  and  glucose  solution  inside  and  outside).  Similar  to  observation 
made  by  (28)  for  spheres,  in  the  single  beam  case,  we  also  notice  that  RO  and  MST 
analysis  matches  exactly  in  the  illuminated  surface  {6  =  7r/2  to  vr  )  for  spheroids 
with  low  eccentricity  due  to  the  lack  of  interference  fringes.  But  in  the  shadowed 
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region  {6  =  n  io  2ti)  the  presence  of  surface  waves  (74)  and  fold  caustics  (75)  (region 
on  the  shadowed  surface  where  rays  from  two  different  angle  of  incidence  strikes) 
causes  wave  interference,  resulting  in  the  deviation  of  RO  from  MST.  The  first  peak 
in  the  RO  stress  profile  arises  from  the  first  order  fold  caustics  (i.e.  on  the  first 
reflection  from  the  shadowed  region)  and  as  is  evident  from  the  figure  it  broadens 
with  eccentricity  due  to  an  enhanced  focusing  effect.  For  higher  eccentric  objects  the 
first  order  fold  caustics  exits  over  a  larger  region,  therefore  the  interference  is  more 
prominent  (extreme  eccentric  case  e  =  0.9).  Fig  5.4(b)  illustrates  the  single  beam 
field  profile  for  a  sphere  (e  =  0)  and  highly  eccentric  prolate  spheroid  (e  =  0.9) 
purporting  the  fact  of  enhanced  focusing  and  interference  effects.  The  dual  beam 
case  (Fig  5.4(c))  is  just  an  extension  of  the  single  beam  stress  profile  and  as  expected 
the  interference  effects  now  become  predominant  throughout  the  surface  of  the  spheroid. 

The  correlation  between  RO  and  MST  analysis  with  refractive  index  variation  and 
size  variation  is  shown  in  Fig  5.5.  The  values  of  refractive  index  contrast  chosen  here 
agree  with  those  which  exist  for  GUVs  (An  =  0.0105),  RBCs  (An  =  0.045)  and 
polystyrene  (An  =  0.245).  As  refractive  index  contrast  is  increased,  the  angle  of 
refraction  increases  and  rays  undergo  larger  deviation  from  the  incidence  direction; 
consequently  the  two  ray  inference  zone  (first  order  fold  caustics)  broadens.  This  is 
evident  in  Fig  5.5(a)  where  the  first  peak  in  the  RO  analysis  becomes  wider  as  we 
increase  the  refractive  index  contrast.  This  trend  leads  to  enhanced  interference  effects 
in  the  shadowed  region  and  thus  causing  RO  to  further  deviate  from  FDTD  analysis. 
However,  in  case  of  size  variation,  the  interference  effects  increase  marginally  with  size 
and  thus  the  variation  between  RO  from  MST  remains  constant. 
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Figure  5.5:  Comparison  of  MST  and  RO  for  spheroids  with  e  =  0.4  (a)  refractive  index  variation 
and  (b)  size  variation. 

However,  in  all  the  eases  studied,  it  is  quite  apparent  that  the  area  under  the  stress 
eurves  for  both  the  RO  and  MST  remains  approximately  equivalent.  Therefore  it  is 
justified  in  using  the  RO  approximation  while  ealculating  the  total  foree  on  the  objeet. 
However,  for  the  stress  profile  calculation,  the  interference  effects  on  the  spheroid  sur¬ 
face  causes  the  RO  and  MST  to  deviate  from  one  another  substantially.  Resultantly, 
for  cases  where  the  point  of  interest  is  surface  stress  (as  in  the  case  for  living  cells  and 
GUVs  trapped  in  DBOT)  MST  provides  a  more  accurate  result  even  if  the  spheroid  size 
is  within  the  RO  limit. 


5.2  Dual-beam  optical  trap  setup 

A  schematic  of  a  dual-beam  optical  trap  is  shown  in  Fig  5.6.  Light  is  emitted  from  two 
optical  fibers  (OF)  creating  a  counter-propagating  Gaussian  beam  trap.  The  optical 
fiber  used  is  HI780  (Corning,  USA)  with  ~  2.4/im  mode-field  radius.  The  beam  size 
at  the  center  of  the  flow  channel  (wq  =  11.5/xm)  was  calculated  from  the  size  of  the 
optical  fiber  mode  using  the  ABCD  matrix  method.  Radiation  pressure  pulls  the  GUV 
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Figure  5.6:  (a)  Schematic  of  the  stretching  of  a  GUV  using  DBOT,  (b)  Optical  setup  that  incor¬ 
porates  a  silicon  chip  for  fiber-to-capillary  alignment  and  microliuidic  adapters  to  couple  the 
flow  channel. 


into  the  center  of  the  two  beams,  creating  a  stable  trap  (68;  69).  In  this  position,  the 
surface  stress  (70;  72)  can  deform  the  GUV  along  the  beam  axis,  as  shown. 


The  micro-fluidics  channel  was  created  with  a  square  capillary  of  100  /xm  inner 
diameter  (Vitrocom,  USA)  and  a  100  /xm  wall  thickness.  The  optical  fibers  were 
aligned  perpendicularly  to  the  micro-fluidics  channel  at  a  distance  of  300  /xm.  The  axis 
of  each  optical  fiber  was  aligned  so  that  the  trapping  region  was  at  the  center  of  the 
flow  channel;  thus  trapped  GUVs  were  not  in  physical  contact  with  capillary  walls. 
The  capillary-to-fiber  union  was  coated  with  index  matching  liquid,  and  the  capillary  is 
coupled  to  a  peristaltic  pump  (Instech,  USA)  using  microfluidic  adapters  (Upchurch, 
USA). 


A  wavelength  of  808  nm  was  used  and  a  power  output  of  250  mW  was  measured 
at  each  fiber.  The  wavelength  was  selected  to  minimize  water  absorption  (76)  in  the 
infrared  and  thereby  avoid  heating  effects  (64;  65).  For  a  wavelength  of  808  nm, 
we  have  simulated  the  heating  of  our  flow  channel  in  COMSOL  and  found  that  the 
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temperature  increases  less  than  1  K  for  250  mW  of  laser  power. 

A  microscope  was  used  to  observe  the  experiment  and  a  video  was  recorded  at  61  fps 
by  a  GiGe  camera  with  CCD  image  sensor  (Easier  AG,  Germany).  Both  the  laser  and 
camera  were  controlled  through  a  customized  Labview  program  (National  Instruments). 

5.3  GUV  fabrication 

To  facilitate  trapping,  it  is  necessary  to  create  a  refractive  index  difference  between  the 
interior  and  exterior  of  the  GUV.  To  achieve  this,  we  fabricated  GUVs  with  an  internal 
sucrose  solution  and  transferred  them  to  an  osmotically-balanced  glucose  solution.  The 
density  difference  between  sucrose  and  glucose  results  in  a  refractive  index  difference. 

GUVs  with  sucrose  inside  were  fabricated  by  electroformation  from  pure 
dioleoylphosphatidylcholine  (DOPC)  using  a  modified  vertion  of  (77).  The  synthetic 
lipid,  l-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine  (POPC),  1,2-dipalmitoyl-sn- 
glycero-3-phosphocholine  (DPPC),  and  cholesterol  were  purchased  from  Avanti  Polar 
Lipids  (Alabaster,  AL).  The  lipids  were  dissolved  in  chloroform  and  deposited  onto  an 
indium-tin  oxide  (ITO)  coated  glass  slide  from  Delta  Technologies  (Loveland,  CO). 
The  electroformation  chamber  was  formed  with  two  ITO  slides,  with  the  conducting 
sides  facing  inwards  and  separated  by  a  2.5  mm  thick  silicone  spacer.  A  500  mM 
sucrose  solution  in  4  mM  HEPES  titrated  to  pH  7.0  with  sodium  hydroxide  was  added 
to  the  overnight  vacuum-dried  lipid  film  on  the  slides  giving  a  final  lipid  concentration 
of  0.25  mg/mE.  A  2.65  V  AC  electric  field,  generated  at  100  Hz  by  a  function  generator, 
was  applied  to  the  chamber.  The  electroformation  temperatures  were  chosen  so  that 
the  lipids  remain  in  the  liquid  phase.  Transfer  to  a  glucose  solution  was  achieved  by 
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performing  three  two-fold  serial  dilutions  on  the  GUV  suspension  using  an  identical 
buffer  solution,  but  with  the  sucrose  swapped  for  an  equimolar  quantity  of  glucose. 

Refractive  indices  were  measured  for  the  pure  sucrose  and  glucose  solutions  using  a 
refractometer  (PAL-RI  Refractometer  -  Atago,  Bellevue,  WA).  Values  were  obtained  for 
the  pure  glucose  solution  (RI  =  1.3455  ±  0.0003)  and  sucrose  solution  (GUV  interior 
RI  =  1.3575  ±  0.0003),  and  the  dilution  used  in  the  experiment  yielded  a  final  refractive 
index  difference  of  0.0105  ±  0.0003,  slightly  less  than  the  difference  between  the  two 
pure  solutions.  The  prepared  mixture  was  then  left  to  sit  for  ten  minutes  in  order  to 
allow  GUVs  to  sediment  to  the  bottom  of  the  tube.  Following  this  time,  half  of  the  total 
liquid  volume  was  carefully  removed  at  the  surface.  This  suspension  was  then  directly 
pumped  into  our  microfluidic  device. 

5.4  Trapping  and  stretching  of  GUVs 

GUVs  were  captured  using  minimal  power  (total  power  of  100  mW;  50  mW  from  each 
fiber)  while  ensuring  the  flow  was  stopped  and  the  GUV  was  resting  at  the  bottom  of 
the  channel.  As  a  GUV  is  pulled  up  toward  the  optical  axis,  its  initial  circular  shape 
is  drawn  into  a  slightly  prolate  elliptical  shape  due  to  the  stress  profile  created  by  the 
optical  stress.  Axial  deformation  of  a  GUV  made  of  POPC  lipid  is  shown  in  Fig  5.7(a). 
The  total  power  was  then  increased  to  500  mW,  and  the  deformed  shape  at  maximum 
power  is  shown  in  Fig  5.7(b).  The  major  axis  is  increased  from  d  =  1 1.53  ±  0.05  /rm  to 
d  =  11.94  ±  0.05  /xm  along  the  beam  axis,  while  the  minor  axis  decreased  from  10.29 
±  0.05  fxm  to  10.05  ±  0.05  fxm.  Fig  5.7(c)  shows  the  contours  at  both  minimum  and 
maximum  power.  Stretching  of  the  GUV  along  the  beam  axis  can  be  clearly  observed. 
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Figure  5.7:  (a)  POPC-GUV  optically  trap  at  low  power,  (b)  POPC-GUV  optically  trap  at  high 
power,  (c)  A  plot  of  the  contours  fitted  to  both  stretching  powers  (blue  =  low  power  /  low  tension, 
green  =  high  power  /  high  tension).  The  scale  bar  is  10  /xm. 

5.4.1  Calibration  of  optical  forces 

The  optical  power  calibration  of  our  DBOT  setup  was  done  using  shooting  experiments. 
We  first  trapped  a  polystyrene  bead  (n  =  1.615)  in  water  with  equal  power  from  both 
beams  and  pushed  it  towards  the  one  beam  by  increasing  the  power  of  the  opposite 
beam.  We  then  blocked  on  the  opposite  beam  by  reducing  the  current  supplied  to 
zero.  The  polystyrene  bead  accelerated  towards  the  opposite  beam  and  the  distance 
vs.  time  information  were  recorded.  The  velocity  as  a  function  of  distance  was  then 
calculated  and  used  in  drag  force  equation,  F^rag  =  QTrripv  ,  where  p  is  the  viscosity 
of  the  water  (0.0009  Pa-s),  p  is  the  radius  of  polystyrene  bead  (4.97  pm)  and  v  is  the 
velocity  of  the  particle.  The  following  correction  factor  was  multiplied  to  the  drag  force. 


c  = 


A  (£)  I  1  (pV 

16  \  b)  ^  8  \  b) 


JA  (P)^ 

256  \  bJ 


where  b  is  the  distance  to  the  closest  wall  and 


2Z/  +  6/3 
_  A  3u  +  b(3 

I3~^  is  the  slip  coefficient. 


(5.4) 
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Figure  5.8:  The  drag  forces  in  no-slip  and  perfect-slip  cases  are  plotted  against  the  calculated 
optical  forces  (solid  lines)  from  one  beam  for  two  power  levels,  50mW  and  lOOmW.  50  /um 
corresponds  to  the  center  of  the  channel. 

The  first  term  is  assoeiated  with  the  proximity  of  the  walls  while  the  seeond  term 
takes  into  aeeount  the  slip-flow  boundary  eonditions  on  the  surfaee  of  the  partiele  [36]. 
The  experiment  was  repeated  with  two  different  optieal  power  levels,  50  mW  and  100 
mW.  The  eorreeted  drag  forees,  assuming  no-slip  eonditions  (/3  — )■  cxd  ),  are  the  upper 
boundary  of  the  foree  exerted  on  the  bead.  In  the  ease  of  perfeet  slip  ( /3  — )■  0),  the 
eorreeted  drag  foree  is  the  minimum  foree  expeeted  on  the  bead  and  therefore  our  lower 
boundary.  In  Fig  5.8  these  boundaries  were  plotted  against  the  theoretieal  optieal  forees. 
Our  theoretieal  predietion  of  optieal  foree  is  eonsistent  within  these  boundaries  while  the 
upper  boundary  matehes  the  theoretieal  optieal  foree  within  ~15%  at  the  eenter  of  the 
ehannel. 

5.4.2  Temporal  response  to  applied  stress 

We  measured  the  response  of  a  GUV  to  a  step  inerease  in  applied  stress.  The  total  laser 
power  was  inereased  from  100  mW  to  500  mW,  as  shown  in  Fig  5.9(a)  (blue  line;  right 
axis). 
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Figure  5.9:  (a)  The  optical  power  (blue  line;  right  axis)  is  suddenly  increased  from  100  mW  to 
500  mW.  The  major  axis  strain  is  shown  by  the  red  dots  (left  axis),  (b)  Micrograph  of  deforming 
GUV.  The  scale  bar  is  10  /rm. 

The  power  was  held  at  its  maximum  value  for  5  seconds  and  then  decreased  to  its 
initial  value.  The  major  axis  strain  is  shown  on  the  left  axis  (red  dots).  The  major  axis 
strain  was  calculated  by  assuming  that  the  shape  of  the  GUV  at  maximum  power  is  a 
prolate  spheroid.  We  take  the  diameter  of  a  sphere  with  the  same  volume  as  the  zero- 
power  value  of  the  major  axis.  The  major  axis  strain  is  the  percent  change  in  major  axis 
compared  to  the  zero-power  value.  From  Fig  5.9(a),  it  can  be  seen  that  the  major  axis 
strain  increases  nearly  instantaneously  with  the  step  increase  in  power.  The  initial  strain 
of  8.2  ±  0.4%  increases  by  4.1  ±  0.25%.  Based  on  our  frame  rate  of  61  fps,  we  are  able 
to  capture  2-3  data  points  in  the  transition  region  between  power  levels. 

5.5  Measurement  of  lipid  bilayer  bending  modulus 

In  our  experiment,  we  gradually  ramped  the  laser  power  from  100  mW  to  500  mW  as 
shown  in  Fig  5.10(a).  The  power  was  increased  in  11  steps,  with  a  holding  time  of  1  s 
at  each  step.  As  the  power  increased,  we  observed  axial  deformation  of  the  GUV  and  a 
clear  dampening  of  membrane  fluctuations. 
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Figure  5. 10:  Measurement  of  the  bending  modulus  of  a  pure  POPC  GUV.  Time  axis  is  common 
to  plots  a-c,  dotted  guide  lines  show  time  points  at  which  laser  power  is  increased,  (a)  Laser 
power  as  a  function  of  time,  (b)  2D  contour  plot  showing  the  radius  as  a  function  of  angle  in  the 
image  plane  as  a  function  of  time,  (c)  Percentage  area  strain  as  a  function  of  time,  (d)  Average 
stress  on  the  GUV  as  a  function  of  eccentricity  and  base  radius,  (e)  Average  percentage  area 
strain  for  each  laser  power  plotted  versus  the  scaled  lateral  tension,  k  is  the  fitted  value  of  the 
bending  modulus. 


We  calculate  the  area  strain  of  a  GUV  as  follows.  The  grayscale  image  for  each 
frame  was  processed  using  our  in-house  Matlab  code  to  trace  the  edge  of  the  GUV.  The 
nominal  edge  is  located  at  the  grayscale  center-of-mass  in  an  8-pixel  region  surrounding 
the  darkest  pixel  at  each  point  along  the  GUV  circumference.  The  edge  coordinates 
formed  a  closed  contour  for  each  frame.  The  nominal  GUV  center  was  located  at  the 
center  of  this  contour. 


The  resulting  contour  data  set  for  each  image  frame  consists  of  the  location  of  the 
contour  center  and  a  contour  radius  for  each  of  360  azimuthal  angular  points.  Contour 
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radius  can  be  expressed  as  a  function  of  azimuthal  angle  (Fig  5.10(b)).  We  then  pro¬ 
ceeded  to  expand  the  contour  data  in  Legendre  polynomials.  This  strategy  fits  the  shape 
of  the  contour  as  an  expansion  of  the  equatorial  plane  of  a  base  sphere  with  radius  R. 
This  base  sphere  is  the  assumed  relaxed  shape  of  the  untrapped  GUV.  To  find  R,  the 
GUV  at  minimum  trapping  power  is  assumed  to  be  a  prolate  sphereoid.  The  volume  of 
this  GUV  is  calculated  as  Vo  =  4/37ra^6,  where  a  is  the  semi-minor  axis  and  b  is  the 
semi-major  axis.  R  is  then  the  radius  of  a  sphere  with  the  same  volume  as  this  prolate 
spheroid.  Initial  area  is  calculated  as  the  surface  area  of  this  sphere,  Aq  =  We 

restrict  our  Legendre  expansion  to  the  second  mode  coefficient  (U2),  since  this  can  be 
used  to  represent  the  shape  change  of  the  contour  from  circle  to  increasingly  eccentric 
ellipse  (78;  79).  We  reconstructed  the  contour  using  the  equation, 

r{e)  =  R  +  U2P2{cose)  (5.5) 


where 


P2{cos6)  =  -{3cos^6  —  1)  (5.6) 

and  U2,  using  the  ortho-normality  of  the  Legendre  polynomials,  is  given  by  the  integral: 

M2  =  ^  y  P2{cos9)[r{cos9)  —  R]d9  (5.7) 

After  reconstruction  the  contour,  the  top  half  the  contour  [O-vr]  was  numerically  inte¬ 
grated  to  generate  a  solid  of  rotation  about  its  major  axis  corresponding  to  the  surface 
area  of  the  GUV 

PTZ 

A  =  2tt  j  rsin9 

Jo 


{d9y 


(5.8) 
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Area  strain  is  calculated  using  {A  —  Aq)/Aq  (Fig  5.10(c)).  As  expected,  we  observe  an 
increase  in  area  strain  with  increasing  laser  power. 


The  bending  modulus  of  the  GUV  membrane  can  be  obtained  by  calculating  area 
strain  as  a  function  of  lateral  tension.  In  the  low-stress  regime  (80) 


^  -  ^0  _  kT 

Ao  Sttkb  V  o-Q 


(5.9) 


where  k  is  Boltzmanns  constant,  T  is  temperature,  ah  is  the  lateral  tension  on  the 
membrane,  and  do  is  the  intrinsic  lateral  tension. 


In  order  to  determine  the  lateral  tension  on  the  membrane  at  each  power  level 
(81),  it  is  necessary  to  calculate  the  surface  stress  on  the  GUV.  Ray  optics  approaches 
have  previously  been  used  to  calculate  the  force  on  spherical  (82)  and  spheroidal  (83) 
objects.  We  assume  a  spheroidal  shape  for  the  GUV,  and  calculate  the  total  force  on  the 
front  and  back  surfaces.  For  each  power  level,  we  calculated  the  force  on  a  spheroid 
with  major  and  minor  axes  equal  to  the  average  values  over  all  image  frames.  We 
included  the  effect  of  multiple  reflections,  up  to  5  bounces.  For  each  incident  ray  and 
each  bounce,  we  determine  whether  the  bounce  occurs  on  the  front  or  back  surface  and 
store  the  vector  force.  For  each  surface,  we  then  add  the  force  contributions  vectorially 
to  determine  the  total  force  on  the  surface.  The  stress  is  calculated  by  dividing  the  total 
force  by  the  surface  area. 

The  calculated  average  stress  is  shown  in  Fig  5.10(d).  We  present  the  results  as 
a  function  of  eccentricity  (e)  and  base  radius  (R).  The  optical  power  from  each  beam 
was  taken  to  be  250  mW  and  the  refractive  index  difference  (An)  to  be  0.0105.  While 
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the  average  stress  decreases  slightly  with  eccentricity,  it  substantially  decreases  with 
increasing  base  radius  (R),  in  the  8-14  /rm  range. 

The  average  stress  is  then  translated  to  a  lateral  tension  ah  (81;  78).  After  each 
increase  in  optical  power,  the  GUV  reaches  an  equilibrium  shape  with  a  homogeneous 
surface  tension  in  the  membrane.  The  pressure  difference  created  inside  (Ap)  is  equal 
on  the  equator  and  the  poles  of  the  ellipsoid.  If  we  assume  that  there  is  no  applied  stress 
at  the  poles,  given  the  ray  optics  approximation,  we  can  use  the  Young-Laplace  equation 
to  write 


Ap  (Ci  T  C2)equ(^h  (Ci  T  C2)pol^h  (5.10) 

where  ci  and  c-2  are  the  principal  curvature  values  at  both  the  equator  and  the  poles,  and 
Tequ  is  the  magnitude  of  stress  exerted  on  the  surface  at  the  equator.  We  calculated  the 
curvatures  from  the  extracted  contours.  The  optical  stress  Tf.qu  was  calculated  using  ray 
optics.  As  mentioned  before,  the  surface  stress  is  zero  at  the  poles  and  the  symmetrical 
discontinuities  or  spikes  in  the  stress  profile  around  the  poles  are  a  result  of  multiple 
internal  reflections.  In  our  calculation,  we  included  the  effects  of  multiple  bounces 
(up  to  5)  of  each  incident  ray  inside  the  GUV.  For  the  refractive  indices  of  the  inte¬ 
rior  and  exterior  of  the  GUV  used  in  our  experiment,  the  reflection  is  small  (R=0.002%). 

We  calculate  an  initial  tension  (do)  of  5.76  ±  0.25  x  10“®  mN/m  and  plot  area 
strain  as  a  function  of  the  log  of  scaled  lateral  tension  in  Fig  5.10(e).  The  error  bars  on 
both  axes  are  equal  to  the  standard  deviation  of  the  corresponding  quantity,  taken  over 
all  images  recorded  at  a  fixed  laser  power.  The  slope  is  proportional  to  the  bending 
modulus,  which  is  found  to  be  7.95  ±  0.45  kT.  The  log-linear  relationship  indicates  that 
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we  are  in  the  low  stress  regime  and  that  area  expansion  of  the  membrane  comes  from 
damping  bending  fluctuations,  as  opposed  to  direct  stretching  (i.e.  area  dilation)  of  the 
membrane,  as  observed  at  higher  stresses  (84). 

We  note  that  the  experimental  data  shown  in  Figs  5.10(a),  5.10(b),  5.10(c),  and 
5.10(e)  is  obtained  from  a  single  GUV.  Moreover,  we  note  that  since  the  stress  is  not 
uniform  over  the  GUV  surface,  a  more  sophisticated  model  of  vesicle  deformation 
would  include  the  effects  of  stress  non- uniformity  on  final  shape.  This  is  an  interesting 
area  for  further  research. 


5.6  POPC,  POPC-col  and  DPPC  bending  modulus 

Under  experimental  conditions  (at  room  temperature),  both  pure  POPC  and  POPC  with 
20  mol%  cholesterol  are  in  the  liquid  phase  whereas  DPPC  with  20  mol%  cholesterol  is 
in  the  gel  phase  (85). 

The  bending  modulus  scales  with  membrane  thickness  squared  multiplied  by  the 
area  compressibility  modulus  (k  ~  h'^Ka  )(86).  Gel-phase  bilayers  are  thicker  (87),  and 
the  bending  modulus  of  a  lipid  bilayer  in  the  gel  phase  is  consequentially  higher,  than 
that  of  the  same  membrane  in  the  liquid  phase  (88;  89). 

We  first  demonstrated  that  optical  stretching  is  able  to  measure  the  difference  in 
bending  modulus  of  lipid  membranes  with  different  phase  states.  Fig  5.11  presents  the 
dependence  of  relative  area  strain  on  applied  tension  for  two  GUV  compositions:  POPC 
and  DPPC  with  20  mol%  cholesterol.  The  initial  tensions  are  at  7.53  ±  0.73  x  10“^ 
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Figure  5.11:  Apparent  area  strain  versus  In(tension)  plots  for  POPC  and  DPPC  GUVs  with 
20%  cholesterol.  The  POPC  lipid  membrane  is  more  flexible  than  DPPC  lipid  membrane.  The 
error  bars  represent  one  standard  deviation  of  area  strain  and  lateral  tension.  The  errors  in  bend¬ 
ing  modulus  values  are  calculated  from  the  standard  errors  of  the  slope  using  linear  regression 
analysis. 

mN/m  and  1.61  ±  0.19  X  10  ^  mN/m,  and  the  calculated  bending  modulus  values  are 
9.16  ±  0.72  kT  and  24.89  ±  2.47  kT  for  the  POPC  and  DPPC  membrane  respectively. 

To  obtain  population  statistics,  we  measured  the  bending  moduli  of  GUVs  fabricated 
from  pure  POPC  (n  =  56),  POPC  with  20  mol%  cholesterol  (n  =  40),  and  DPPC  with  20 
mol%  cholesterol  (n  =  39)  respectively.  We  made  fit  normal  distributions  to  the  data  to 
obtain  the  following  bending  modulus  means  and  standard  deviations:  8.13  ±  2.06  kT 
for  POPC  vesicles,  8.50  ±  1.83  kT  for  POPC  with  20%Chol,  and  27.24  ±  7.69  kT  for 
DPPC  with  20%Chol.  Histograms  of  the  population  bending  modulus  data  are  given  in 
Fig  5.12. 

We  confirmed  the  normality  of  the  populations  collected  using  the  Kolmogorov- 
Smimov  test.  We  then  performed  the  students  t-test  to  compare  the  two  POPC 
populations  (POPC  and  POPC-20%Chol).  The  populations  are  not  significantly  differ¬ 
ent  (p  >  0.05).  This  analysis  confirmed  the  indistinguishability  of  POPC  populations 
with  and  without  cholesterol,  in  contrast  to  other  published  results  (90;  91).  Cholesterol 
is  believed  to  increase  n  by  ordering  the  acyl  chains.  Marsh  did  an  extensive  comparison 
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Bending  Modulus  (kT) 

Figure  5.12:  Histograms  of  the  bending  modulus  distributions  of  and  DPPC-Chol,  POPC-Chol, 
and  pure  POPC  lipid  GUVs. 

of  several  methods  and  eoneluded  that  eholesterol  inereases  mean  bending  moduli  (92). 
However  reeent  reports  have  argued  that  the  eontribution  of  eholesterol  to  bending 
modulus  is  not  universal  and  unsaturated  lipids  (partieularly  1,2-dioleoyl-sn-glyeero- 
3-phosphocoline,  DOPC)  are  not  affected  by  the  addition  of  cholesterol  (78;  93).  Our 
results  agree  with  these  recent  reports.  We  were  unable  to  produce  high  yields  of  pure 
DPPC  GUVs. 

There  are  several  factors  that  may  influence  the  width  of  the  distributions  shown 
in  Fig  5.12.  GUVs  made  from  POPC  may  be  photooxidized  when  exposed  to  bright 
light  during  microscopy  (94).  However,  in  our  optical  stretching  experiment,  the 
experimental  duration  for  a  single  GUV  was  around  30s,  and  no  fluorescent  dye 
was  used,  so  we  do  not  expect  significant  photooxidation  to  occur.  The  width  of  the 
DPPC/cholesterol  distribution  may  also  be  a  result  of  compositional  heterogeneity  due 
to  non-uniform  cholesterol  incorporation  during  lipid  film  drying,  as  has  been  shown 
previously  (95). 
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We  analyzed  possible  sources  of  systemic  error  in  this  experiment.  One  potential 
error  source  comes  from  how  we  calculated  vesicle  surface  area,  by  fitting  the  thermal 
undulations.  For  the  data  presented  in  Fig  5.12,  only  the  second  Legendre  coefficient 
is  used  to  calculate  the  area  change.  When  we  include  the  fourth  mode,  as  suggested 
in  by  Gracia  and  coworkers  (78),  the  bending  modulus  values  tend  to  decrease  by 
approximately  5%,  within  the  reported  error.  Including  the  fourth  mode  does  not 
decrease  the  width  of  the  population  histograms,  indicating  that  lack  of  precision  in 
shape  fitting  is  not  a  major  contributor  to  the  width  of  the  observed  distributions. 

Second,  GUVs  are  trapped  at  a  position  slightly  below  the  axis  of  the  optical  trap 
due  to  gravity.  For  example,  consider  a  spheroidal  GUV  with  10  /rm  equivalent  sphere 
volume  and  eccentricity  of  0.4.  The  net  gravitational  force  on  the  GUV  is  proportional 
to  the  density  difference  between  the  solutions  on  the  inside  and  outside  and  equals 
~1.3  pN.  The  trapping  position  can  be  calculated  by  balancing  the  gravitational  force 
with  the  optical  gradient  force.  At  the  minimum  trapping  power  of  100  mW,  the  GUV 
is  trapped  ~2  /rm  below  the  optical  axis.  At  a  power  of  500  mW,  the  offset  decreases  to 
~  0.4  /rm.  The  microscope  is  initially  focused  on  the  equatorial  plane  of  the  GUV.  As 
the  GUV  is  pulled  further  upward  into  the  trap,  the  segment  of  the  vesicle  in  the  focal 
plane  will  shift. 

To  test  whether  this  effect  could  have  a  significant  impact  on  the  results,  we 
calculated  the  volume  of  the  GUV  at  every  frame  by  assuming  the  captured  contour 
corresponded  to  the  equatorial  plane  of  a  spheroidal  vesicle  and  measuring  the  length  of 
the  major  and  minor  axes.  If  the  upward  shift  of  the  GUV  were  significant,  it  would  be 
expected  to  result  in  consistently  smaller  measured  GUV  volumes  as  power  increased. 
This  is  because  the  true  volume  of  the  GUV  (expected  to  remain  constant  through  this 
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experiment)  can  be  calculated  from  the  true  equatorial  plane  contour  measured  when  the 
microscope  is  initially  focused  at  minimum  trapping  power.  Contours  corresponding 
to  other  GUV  cross-sections  would  result  in  a  smaller  calculated  volume.  A  significant 
shift  of  the  vesicle  along  the  optical  axis  would  therefore  result  in  an  apparent  (and 
purely  artifactual)  decrease  in  volume.  In  fact,  in  no  experiment  did  the  average 
calculated  volume  change  exceed  1%,  and  there  was  no  relationship  between  volume 
and  power. 

Finally,  we  performed  a  set  of  experiments  to  test  whether  the  GUV  stretching 
system  exhibits  hysteresis.  The  optical  power  was  ramped  up  and  down  between 
minimum  and  maximum  trapping  power.  The  bending  modulus  was  obtained  both  as 
the  vesicle  dilated  and  contracted.  Fig  5.13  shows  two  examples  of  time-dependent 
area  strain  and  extracted  area  strain  vs.  surface  tension,  for  pure  POPC  (Fig  5.13(a)) 
and  DPPC-20%Chol  (Fig  5.13(b)).  For  this  POPC  vesicle,  the  bending  modulus  values 
obtained  from  the  two  slopes  gave  6.52  ±  0.32  kT  for  ramp-up  and  6.58  ±  0.35  kT 
for  ramp-down.  For  this  DPPC-20%Chol  vesicle,  we  obtained  27.47  ±  3.29  kT  for 
ramp-up  and  28.04  ±  2.25  kT  for  ramp-up  and  ramp-down.  In  all  cases  tested,  the 
change  in  measured  bending  modulus  for  the  two  different  deformation  directions  is 
within  error. 

In  conclusion,  we  have  used  optical  stretching  in  a  microfluidic  DBOT  system  to 
investigate  the  elastic  bending  modulus  of  populations  of  POPC,  POPC  with  20%  Choi, 
and  DPPC  with  20%  Choi  GUVs.  The  difference  in  bending  modulus  between  liquid 
and  gel-phase  GUVs  is  unambiguous;  gel-phase  GUVs  are  significantly  stiff er.  The 
bending  elastic  modulus  of  the  membrane  composed  of  the  monounsaturated  phospho¬ 
lipid  (POPC)  is  insensitive  to  Choi  at  concentrations  up  to  20%.  This  technique  is  a 


72 


84 


»»  yiJ! 

::f !' 
$  8®^ 
io  ■=*  < 

OJ 

"••••U*  Ij 

M«nipOo«n 

X 

0  S  10  IS  20 
T«ne  (s) 

-7«  -72  -«8  -« 
In(o-mNAn) 

"  5" 

Oil  |os 

Rwieue  ^ 

r  ■ 

0  S  to  15  20 
Time  (8) 

-7 

S  -72  -«4  -4 

Figure  5.13:  Laser  power  ramp-up  and  ramp-down  experiments  on  two  different  GUVs,  (a) 
POPC  and  (b)  DPPC-20%Chol. 


simple  and  powerful  approaeh  to  the  measurement  of  membrane  bending  properties.  It 
has  the  potential  to  be  widely  deployed  in  efforts  to  understand  relationships  between 
membrane  composition  and  membrane  mechanics. 
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